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Taylor Diffusion in Tubular Reactors’ 


JOSEPH COSTE*, DALE RUDD® and 
NEAL R. AMUNDSON* 


In order to estimate the effect of velocity profile 
and radial diffusion on conversion in tubular reactors 
the assumption of Taylor diffusion is adopted. This 
leads to equations which are simple, but because of 
the numerical instability of these equations they are 
approximated to by a cell model in order to obtain 
numerical results. It is shown that the cell model 
can be a very good approximation to the one dimen- 
sional diffusion equation. If the Taylor diffusion 
model is valid for systems involving chemical reaction 
and heat transfer it may be concluded that while 
there seem to be no gross effects in a tubular turbu- 
lent reactor there are differences in detail, such as 
maximum temperature at the hot spot, which may be 
important. 


[" order to solve the problem of a tubular reactor to determine 
the conversion and heat transfer characteristics it is necessary 
to consider a rather complicated system of partial differential 
equations. This system is complicated because the velocity 
profile is not flat and because there are gradients of concentration, 
temperature, and pressure. To solve ‘this problem in a general 
way analytically is, of course, out of the question and the 
numerical solution is tedious and requires a very large computer. 
At the same time there is the larger question whether from a 
practical point of view the effects mentioned above are important 
enough to justify the effort needed to solve the problem. It 
would be interesting, however, to obtain some estimate of the 
magnitude of some of these effects, and it is the purpose of this 
paper to exploit a technique which might give a partial insight 
into the problem. 

G. |. Taylor“:» and later Sjenitzer“ and Tichacek, 
Barkelew, and Baron“) showed that the effect of nonflat 
velocity profile and radial mixing under the assumption of the 
Reynolds’ analogy could be combined into a single parameter, 
a quasi longitudinal mixing coefficient which appears in the 
differential equations as a diffusivity. This has the advantage 
that the partial differential equation is replaced by an ordinary 
differential equation whose solution is more tractable. The 
above mentioned authors considered the problem in which only 
two modes of transport were important, that by the non- flat 
profile which causes gradients of concentration and the resultant 
transport by radial diffusion caused by the gradients set up by 
the profile. Other modes of eee were considered to be 
negligible compared with these. Taylor showed that in laminar 
flow the proper quasi longitudinal diffusion coefficient was 


a? V2 
48D 


where D is the molecular diffusivity, while for turbulent flow 
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D* = 104eV ./—— 
pV? 


where Vr/pV? is a dimensionless friction factor. 


The basic assumption of this paper is that the quasi mixing 
coefficient of Taylor may be used in the tubular reactor in which 
chemical reaction, heat transfer, and diffusion are coupled. 
While it is undoubtedly true that the assumptions under which 
Taylor diffusion was derived are not valid in this case, it is 
hoped that an estimate of the effect of profile and radial diffusion 
might be obtained for comparison with the lumped constant 
reactor. 


General Problem 


For the purposes of this paper an incompressible fluid with 
a single irreversible first order exothermic reaction will be 
considered. If ¢ is the concentration of reactant then the mass 
conservation equation with Taylor diffusion is 


~?e _dc —-— 
D dx? — J a = cpe FES acetal (1) 


where V is the average velocity of the turbulent flow field. 
The equation of energy conservation with the usual simplifying 
assumptions 1s 


E 
LT iT - — 
pe. FS +R = The eee RP. 2. 
dx? dx 
where the Reynolds analogy is assumed to hold. With these 
equations must be associated some inlet conditions. For our 
purposes the following will be used, 
dc 
Vein = Ve — D* —; x = 0 
dx 
ie oe dT 
VI; = VT —- »* -¢=@ 
dx 
It is convenient to make the change of variable v = x/V to 
give 
D* @ d z 
“C ac = 
_ =cpe RT.... ~ 63) 
V2 dv? dv — ( 
a) - E 
D*@T dT ae oo aoa 
——+K(T, —T) = —cpe RTQ . (4) 
V2 dv? dv 
a D> dc (5) 
Cin = ae ae 
_ _D* aT 
Tin =T-— V2 . . .(6) 
2 dv 


Another condition 
It will be 


where the variable v is a residence time. 
must be stated to completely specify the solution. 


assumed that the reactor is long enough that a finiteness condition 
may be imposed on all solutions. 
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Figure 1—Sequence of forward integration attempts of the 
diffusion equations : temperature. 


Equations (3), (4), (5), and (6) will be solved by the march- 
ing technique, that is, by starting at the reactor entrance and 
proceeding stepwise down the reactor. In order to initiate the 
calculations the values of dc/dv and dT/dv must be known at 
the entrance. These are unknown but a good approximation 
may be obtained from the solution of the same problem treated 
as a lumped constant system. Integrations were then performed 
using the Gill modification of the Runge-Kutta technique. It 
was found that there were no values of the gradients which 
would produce finite solutions. A change of a single unit in the 
sixth significant figure in the inlet gradient makes the solution 
change from large positive to large negative values as shown 
in Figures 1 and 2. Curves 1 thru 7 on these figures show the 
effect of progressively better estimates of the inlet gradients. 
It should be mentioned that the effect of increment size has 
presumably been removed from these calculations. 


The phenomenon described above is not an unusual one and 
is discussed by Fox“. For linear equations the same phenomenon 
may exist since, in effect, one is trying to guess the inlet gradients 
which will give certain values to the arbitrary constants which 
appear in the general solution. For example, if one term of the 
solution becomes very large for large v then its coefficient must 
be exactly zero. In a computer with finite precision this may 
not be possible. 

The possibility of integrating the equations in the reverse 
direction exists but in this case four values must be estimated, 
c, T, dc/dv, and dT/dv at some point. This does not seem 
feasible. 


Finite Cell Model 


One may imagine a tubular reactor divided into cells, each 
of which resembles a stirred pot, of length Ar. It is clear that 
in the limit as Ax approaches zero the series of stirred pots 
approach a tubular reactor with lumped parameters. On the other 
hand for a non-zero Ax the mixing in each of the cells produces 
an effect like that of axial mixing save for the fact there is no 
back mixing. The mass and energy conservation equations 
may be written as 


E 
ViTe-as — Ts) = —cspe RT, Q Ax + K(Ts — Tw) Ax 
E 
Vics —- az — Cx) = Crepe RT, Ax 


where the subscripts denote the cell in which the evaluation is 
to be made. Now if one desires the cell model and continuous 
model to agree the best in some sense one must have a criterion 
for the choice of Ar. One way to proceed is as follows. Let 
Tz be expanded into Taylor series about 7, and choose Ax 
so that the two models agree through second order derivatives. 


150 


K=0,00333 | /sec 
OPV%2.5 sec 


p*6.566 x10" i/sec 


E/R=11,323 17° 


0.02100} = 7300 em® KYgm-mole Ty #3375 °K 







° 
8 
8 
8 


0.01900;—- 


CONCENTRATION gm. moles /cm> 


0.01800}-- 





0171 
ee i0 20 30 40 
TIME v sec 


Figure 2—Sequence of forward integration attempts of the 
diffusion equations : concentrations. 


Thus, 
dT Ax? @T 
Ve -As = = Ax eee 
. de + 2. det 
dc Ax? d?c 
tome te Ae ee 


The discrete model then becomes 


E 
dT @T Ax | ~ 





Y) ee eee wee \T, = 
dx * ae 2 Pe Re Rt Aire ae 
E 
Vv dc + dc Ax| ~ ~ on 
dx dx? 2] ail ; 


Comparison of these two equations with Equations (3) and (4) 
reveals that the equations will be identical through second 
derivative terms if 


VAx 
_ 
or 
V2Av 
—* 
or 


(Ax)? = 2D*Av 
which bears a resemblance to the Einstein relation. 
If we now think of nAx = x, then 
E 


Ax -sz Ax 
1g Rim > pe RTreQ + > MTs ~ Te)i..c... AB 
Ax -=sm 
Ca—1 = Cr + = pe RT ¢, 


and the solution of these equations may be started by specifying 
T, = Tin and ¢, = Cin. For more complicated kinetic relations 
Equations (7) and (8) would require some iterative approximate 
solution but for the case at hand Equation (8) may be solved 
for ¢,, the latter substituted in Equation (7) with the solution 
for T,, following. 
‘Two cases were computed and these are shown on Figures 
3 and 4+. The parameters for the system are given on the 
figures and the two graphs are for cases 
Ax 
2V 


Ax 
= 5 sec. and — = 10 sec. 
2V 


where these would hold for a typical turbulent flow. The 
small circles represent the values of concentration and temper- 
ature versus residence time v. The curves which appear to 
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Figure 3—Mixing cell, backward integration and plug flow 
profiles. D*/V* = 5. 


pass through these points are mot smoothed curves. It is now 
possible to obtain estimates of the concentration, temperature, 
and their gradients for a large v and so integrate the differential 
equations, Equations (3) and (4), backw ards. This was done 
by taking these values at v = 600 and the equations were 
integrated. The curves which appear to pass through the small 
circles were obtained and, further, the inlet conditions, Equations 
(5) and (6), were satisfied within a distance of one cell, which 
is all that one could expect. It is a remarkable fact that the 
mixing cell model and continuous model agree with such pre- 
cision. Using the calculated values of the inlet gradient and 
concentration from the backward integration of the differential 
equations the computation was reversed to determine if the 
same profile could be generated. The instability manifested 
itself again and curves 5, 6, 9 and 10 were generated on Figures 
3 and 4. 


Also shown on Figures 3 and 4 are curves 1 and 2 for plug 
flow with lumped constants. While the profiles differ the main 
difference lies in the maximum temperature attained at the hot 
spot and the fact that it takes longer to produce the same 
conversion in a diffusion reactor than it does in a plug flow 
reactor. 

The calculations in this paper were performed for a typical 
turbulent flow case and the results on Figures 3 and 4 show the 
effect of the turbulent diffusivity. It has been pointed out 
(ref. (4)) that the Taylor diffusion coefficient is very sensitive 
to profile and hence it is possible that the temperature level in 
the reactor is sensitive to those parameters which affect the 
profile. Some additional calculations were made to determine 
the sensitivity of the backward —— of the differential 
equations to the assumed value at v = 600. A change of 1 4% 
in the concentration at v = 600 piniosed a change of 15% in 
the initial values at 7 These profiles are shown in Fi igare 5. 
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Figure 5—Sensitivity of the backward integration. Profiles 
with respect to the starting value at v = 600. 


The Canadian Journal of Chemical Engineering, August, 1961 






p* exp (24.90768) 
E/R*10,757 1/°K 
Q=7300 cm °K/gm mole 
K 0.0225 \/sec 
Tw=340 °K 










a i Tin*340 °K 020 = 
A ‘ Cin *0.020 gm mole/em? Ss 
3 I § 
S37 \ o6 > 
- oL PLUG FLOW z 
¥ (2) ¢ PLUG FLOW S 
> ses} (3) T DIFFUSION BACKWARD 12 
3 (4) C DIFFUSION BACKWARD g 
S | (5) T DIFFUSION FORWARD 3 
2 (6) C DIFFUSION FORWARD E 
F 35s © MIXING CELL DATA 008 ; 





3. 
355 


80 160 240 320 400 480 560 640 
TIME v sec 


Figure 4—Mixing cell, backward integration and plug flow 
profiles. D/V? = 10. 


Conclusions 


It is assumed in this paper that the model of Taylor diffusion 
may be used in a tubular reactor in which there is chemical 
reaction and heat transfer. If the assumption is valid then it 
can be concluded that there are some effects, mainly in the 
reduction of the maximum temperature at the hot spot and 
increased length for the same conversion. These effects are 
not unexpected. A secondary conclusion is that a cell mixing 
model may be defined which gives profiles i in very good agree- 
ment with the Taylor diffusion model for the reactor. It was 
also shown that if diffusion is coupled with exothermic chemical 
reaction the equations describing the process may become very 
unstable and the cell model may provide a reasonable method 
for the numerical solution. It is rather interesting that while 

particular value of the cell size allows for axial diffusion 
(without back diffusion) a small value of Ax produces the 
lumped constant tubular solution. Hence a single mathematical 
model solves both problems. 


Nomenclature 
A = cross sectional area of tubular reactor 
a = radius of tubular reactor 
‘ = molar concentration 
Cp = heat capacity 
D = molecular diffusivity 
D* = effective Taylor diffusivity 
E = activation energy for chemical reaction 
h = overall heat transfer coefficient at tube wall 
hP 2h 
kK = = 
Alc, aC» 

n = cell number 
p = frequency factor for reaction 

AH 7 ; 
Q =- , heat of reaction parameter 

Cp 
T = reactor temperature 
T, = temperature in nth cell 
T. = ambient reactor temperature 
v = x/V, residence time 
V = average fluid velocity 
x = reactor length variable 
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Mass Transfer With and Without 


Chemical Reaction’ 


J. MARANGOZIS? and A. I. JOHNSON? 


Rates of dissolution of cast solids in a baffled 
agitated vessel were obtained for the systems: Benzoic 
acid-water, o-salicylic acid-water, benzoic acid-22.5% 
acqueous glycerol and benzoic acid-41% aqueous 
glycerol. A new continuous technique was developed 
for the measurements, and the data were empirically 
correlated by the method of least squares with the 
equation 


° 0.65 


kyl MK 


“a2 
- = 0.402 me 
D 


pD M 
with a mean deviation of + 5.8%. 


The present data agreed qualitatively with the 
predictions of the turbulent boundary layer theory. 


Reaction factors, k,,/k,°, were also measured in the 
agitated vessels for the systems: Benzoic acid-aqueous 
sodium hydroxide solutions, benzoic acid-aqueous 
potassium hydroxide solutions, and o-salicylic acid- 
aqueous sodium hydroxide solutions. The data were 
found to depend strongly on the ratio of the reactant 
concentrations, C,/C,,, in a manner quantitatively 
predicted by the boundary layer theory corrected for 
ion-diffusion effects. 


Reaction factors were also measured in dia- 
phragm-type diffusion cells for the system benzoic 
acid-sodium hydroxide. The results provided quanti- 
tative support for the ion-diffusion theory in a stag- 
nant film. 


r l Xhere are numerous examples of processes involving diffusion 
accompanied by a chemical reaction. 


Theoretically, at least, the subject has been treated in the 
literature almost exhaustively and as far as its mathematical 
complexities would allow. However, there are relatively fewer 
works reported which endeavored to investigate the problem 
experimentally. 


In the present work, the problem of a rapid, heterogeneous 
chemical reaction was studied experimentally to gain an under- 
standing of the rates and the mechanism involved in two different 
apparatus: a baffled agitated vessel operating under fully- 
developed turbulent conditions, and a diaphragm- type diffusion 
cell with the reaction occurring in the sintered disc under truly 
stagnant conditions. 


\Manuscript received March 1; accepted April 29, 1961. 

*Department of Chemical Engineering, University of Toronto, Toronto, Ont. 
Contribution from the Department of Chemical Engineering, University of 
Toronto, Toronto, Ont. Dr. Johnson is now with the Department of 
Chemical Engineering, McMaster University, Hamilton, Ont. 
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Historical Background 
The Diffusion Layer Theory 


In 1904 Nernst“ advanced the concept of a saturated diffu- 
sion layer at a solid-liquid interface during the dissolution 


process. He suggested that the mass transfer coefficient, k,°, 
be defined as 
, D 
XL 


Whitman’, and Lewis and Whitman“) extended the diffu- 
sion layer concept to cover two films at a gas-liquid interface. 
The ‘‘two-film theory’? was applied by Hatta“) to the absorp- 
tion of CO. in OH™ solutions. Fi igure | shows the concentration 
gradients according to Hatta for the general, rapid chemical 
reaction 

A+B-— AB... ; ; 6 cCae 


Sherwood and Pigford derived an equation for the mass 
transfer coefficient, k,, namely, 
Dz Dp Cz 


ky = 1 ——]..... (3) 
: XL . Da Cai 


It has been reported’: that solubilities, C4;, and diffu- 
sivities, D4 and Dg, of electrolytes are influenced by the presence 
of other ions. Therefore, in interpreting the rates according to 
theoretical models care should be exercised to include the 
aforementioned effects. 


In 1931 Roller®°" demonstrated that the postulation of a 
stationary liquid film, at solid-liquid interfaces, was unrealistic. 
In accordance, other theoretical models ought to be developed. 


The Penetration and Surface-Renewal Theories 


In 1935 Higbie“® advanced the “penetration model” 
mass transfer by unsteady-state diffusion. He derived an 
equation for k° which was radically different from Equation (1), 

D 1/2 


ki? =2 r one (4) 
Tr 


Danckwerts"® and Sherwood and Pigford“ derived an expres- 
sion for the reaction factor, k,/k,°, for a rapid, bimolecular 
chemical reaction, namely 

ky 1 


= mm, (5) 
ki? —serf (a/Da)!'? 


where a is an integration constant which is a function of Cp/C4i, 
D, and Dg. 

Danckwerts"® postulated a “surface-renewal” mechanism 
in which the surface is continually replaced at a constant rate, 
s, and transfer takes place by unsteady- state diffusion. He 
derived an equation identical to Equation (5) for this case, 
too. Hanratty) extended mathematically the Higbie and 


“ 
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Figure 1—Sketch of concentration gradients according to 
Hatta theory of gas absorption with rapid chemical reaction 
in the liquid phase. 


Danckwerts models, and compared his predictions with the 
experimental data of Lin, Moulton and Putnam». 


The Development of the Boundary Layer Theory 


According to Prandtl“, for turbulent flows in pipes the 
fluid is thought to be divided into three layers: the laminar 
sublayer near the wall; the buffer region next to that; and the 
oy peragge core. Von Karman“? obtained equations relating 

to y*. In turbulent liquid streams, however, Lin, Moulton 
and Putnam“) assumed that there is some eddy motion in the 
sublayer which decreases rapidly to zero at the wall. The 
analogy then of mass and momentum transfer leads to the 
equation 


kr 
— bp = f/2 ; (6) 


where # is the mean fluid velocity, f is the skin friction factor 
and ¢p is a complex function of the Schmidt number“®, 


Sherwood and Ryan“® used Deissler’s analogy®® to derive 
the following equations for k,° and k,/k.° for a rapid reaction 


ashy si dy* $( Sc) 7 
—{- = —_————— = rf, Oe (7) 
Re’ \2 , 1/Se + €/v 7 


and 
kr _ Ce — P(H*, p a a) 
ke Cai lb(%0* , S¢a) — (yr ; Sca)) 
1/2 
where yt = 4 ie (9) 
v\p 


In the special case of yg* < 1, Equation (8) was found“ 


to reduce to the simple form 


ba (Sea\"" , (Sca\" Co 
ki? Sta Scp Cai 


It is to be noted that Equation (10) was also derived by 
Potter, and by Friedlander and Litt®” for /aminar boundary 
flows. 


(10) 


lon Diffusion in Interphase Mass Transfer 


worked out the theory for the 
The theoretical 
rate equations were found to agree satisfactorily with experi- 
ments of diffusion of ions with, through or against other ions“), 
Sherwood and Wei?) combined Equation (3) with the principles 
of the V inograd-McBain theory“? to obtain the rate of mass 
transfer in the presence of ionic ‘effects. A test of these calcula- 
tions with experimental data is made in this paper. The 


Vinograd and McBain“ 
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Figure 2—Agitated vessel and accessories. 


1. glass cylinder 6. tie rods 

2. brass plate 7. cast solid 

3. stainless steel plate 8. baffle 

4. impeller 9. liquid distributer 
5. gasket 10. weir hole 


significance of ion diffusion in boundary layer calculations is 
also demonstrated. 


Experimental Part 


(a) Mass transfer in the agitated vessel. 


The apparatus developed by Johnson and Huang) was 
used. It consisted of an agitated vessel the glass cylinder of 
which (6-in. I.D. and 12-in. high) was tightened between two 
plates as shown in Figure 2. A continuous operation was 
adopted mainly because constant conditions would prevail at 
steady state. Also, the calculations were simpler, easier and 
probably more reliable than if an unsteady-state operation had 
been adopted. 

\ shortcoming of the continuous technique was the difficulty 
encountered in the establishment of steady state as discussed 
below. 

The bottom plate of the agitated vessel was of type 304 
stainless steel in which a 14-in.-deep circular ring depression 
was made. An organic acid was cast in the ring by the method 
described by Johnson and Huang®*. The surface area of the 
cast material was 22.8 cm.2. The impeller was a 3-in.-diameter 
flat blade turbine with six blades made of stainless steel. It 
was driven by a Zero-Max variable motor, model M14, type 

G115-CCW. The speed of rotation was checked with a 
Strobotac and it was found to remain constant within less than 
+1% of the set value over a four-hour period. Four vertical 
baffles mounted on two stainless steel rings were used as shown 
in Figure 2. A %¢-in.-O.D., 7-in.-long stainless steel tube was 
bolted on the bottom plate to serve as an overflow weir. 
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Figure 3—Schematie diagram of mass-transfer apparatus. 
1. solvent reservoirs 5. agitated vessel 

2. constant-head vessel 6. solvent-heating coil 

3. solvent flow meter 7. thermostated bath 

4. zero-max motor 8. volumetric cylinder 


A 14-in.-O.D. stainless steel tube, shaped in the form of a 
ring, sitting on the top frame of the bafHles, was employed to 
distiibute the liquid fed into the tank through several small 
holes staggered along the circumference of the tube. 


The agitated vessel was placed in a constant temperature 
bath held at 25 + 0.1°C. during an experiment. A schematic 
diagram of the assembled mass transfer apparatus is shown in 
Figure 3. The systems investigated were: Benzoic acid-water, 
o-salicylic acid-water, benzoic acid-22.5% aqueous glycerol, 
benzoic acid-41% aqueous glycerol, benzoic acid- aqueous sodium 
hydroxide solutions, benzoic acid- aqueous potassium hydroxide 
solutions, and o-salicylic acid-aqueous hydroxide 
solutions. 


sodium 


Three liters of solvent were brought to 25°C. and added to 
the tank to a level of 7-in. It was soon found that steady-state 
concentrations would be reached only after a considerable length 
of time had elapsed, which was too long to have a smooth 
surface maintained throughout. ‘To overcome this difficulty, 
the solvent added to the tank in the case of the physical dissolu 
tion experiments was a solution of the organic acid at approxi 
mately the anticipated steady-state concentration. For the 
chemical dissolution experiments, the solvent added initially 
was a solution of caustic at about the anticipated concentration 
Considerable trial and error effort was devoted to establishing 
the steady-state value; this should be considered as a shortcoming 
of the continuous operation method 


During a run, solvent from the reservoirs, 


1, lowed through 
the constant head vessel, 


2, the flowmeter, 3, the temperature 
adjusting coil, 6, into the tank, 5, being uniformly distributed 
by the liquid distributor and running slowly along the wall of 
the tank. Samples of the ctiluent liquid were taken regularly 
for analysis by titration. The duration of a run was generally 
90 minutes; longer runs were unreliable due to surface irregu 
larities developed. Dissolution was allowed to proceed to the 
extent of 1 mm. in depth, only if a smooth, uniform surface 
was maintained. If the cast acid showed a rough surface at the 
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Figure 4—Schematic diagram of diaphragm diffusion cell, 


1. diffusion cell 5. jacket 


2. diaphragm 6. overflow tube 
7. rubber joint 


8. rotating shaft 


3. stirrer 
4. magnet 


(T)—Temperature measurement. 


end of a run, the results were discarded. These specifications 
restricted the use of caustic solutions to concentrations below 
0.13 N., which corresponded to a value of the reaction factor 
of about 14 as will be discussed later. 


Viscosities and densities of the glycerol solutions were 
measured according to standard methods. Diffusivities of the 
reactants in solutions of the reaction product were determined 
with diaphragm-type diffusion cells and were reported else- 
*) A typical diffusion cell used in the diffusivity measure- 
ments is shown in Figure 4. 


where 


(b) Reaction factors in diaphragm cells. 


\ new procedure was adopted in this work for the measure- 
ment of reaction factors in diaphragm cells. Essentially, the 
diffusion cell illustrated in Figure 4 was equipped with a system 
for the continuous circulation of benzoic acid solution as shown 
in Figure 5. The solution was pumped from reservoir 1, by 
the stainless-steel, diaphragm pump 2, at a rate of about 
10 ce./min. The constant-head vessel 3, was used to regulate 
the flow into the diffusion cell 4, at a rate of about 8 cc./min. 
by letting in a continuous stream of solution to maintain a 
constant acid concentration in the top compartment. The 
cell was magnetically stirred as shown and thermostated at 
2540.1°C. 

Sodium hydroxide solutions were contained in the bottom 
compartment and the reaction was allowed to proceed until 
20 ~ 25% of the base had been neutralized. A_ preliminary 
diffusion of about two hours was allowed in order to establish 
a steady-state gradient in the pores of the diaphragm which 
was of “Fine” porosity, obtainable from Corning Glass, N.Y., 
U.S.A. After the preliminary diffusion, the solutions in the 
compartments were carefully replaced by fresh solutions and 
an ordinary run was started which usually lasted for one or 
two days. At the end of the experiment the concentrations of 
the solutions in the two compartments were determined by 
titration and reaction factors were calculated from the results. 
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Figure 5—Schematic diagram of apparatus for reaction in 


diaphragm cells. 
1. solution reservoir 4. diaphragm cell 
2. diaphragm pump 5. overflow receptacle 


3. constant-head vessel 


(T)—Temperature measurement. 


Discussion of Results 

Physical mass transfer coefficient in the agitated vessel. 

Experimental and*calculated results were listed in Table 1*. 
Mass transfer coefficients were plotted vs. the stirrer Reynolds 
number in Figure 6. For the first three systems the data showed 
that the mass transfer coefficient, k,.°, is proportional to approxi- 
mately the 24 power of Re up to a Rey nolds number of 60,000. 
Beyond that value, the data for the first two systems show what 


appears to be a transition to a different variation of k,° with 
Re in the range 75,000 < Re < 110,000. 


In the more viscous solutions, it was not possible to reach 
this range of high Reynolds numbers with the present apparatus 
to study the phenomenon further. However, despite the repro- 
ducibility of the results in that region the authors have certain 
doubts as to the true significance of those results. It was visually 
observed that at high turbulence the liquid in the tank contained 
a substantial number of finely dispersed air bubbles. It was 
probable that air entrapped in the liquid might cause a reduction 
of the dissolution rate, for example, by reducing the convective 
effects of turbulence. 


The data for the high Sc number system, benzoic acid-41% 
aqueous glycerol, showed a more complicated dependence pr 
k° on Re. This fact pointed to the conclusion that the relation 
among k;°, Re and Sc groups may not be a simple one over a 
wide range of the variables. For a comparison of the present 
data with other investigations on solid-liquid mass transfer in 
agitated vessels the reader is referred elsewhere®**®, It would 
be extremely desirable to compare the present data directly 
with the predictions of the turbulent boundary layer Equations 
(6) or (7) previously discussed. Unfortunately, there are no 
data available in the literature for skin friction factors, froin 
fully baffled agitated vessels. Therefore, the best that can be 
done at the present is to correlate the data empirically and 
make qualitative comparisons with the theories. 

The Gilliland-Sherwood form of correlation provides, on 


dimensional analysis basis, a satisfactory representation of mass 
*Tables 1, 2 and 3 of this paper have been deposited as Document No. 
6718 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be 
obtained by citing the Document No. and by remitting $5.00 for photo- 
prints, or $2.25 for 35 mm. microfilm, Advance payment is required. 
Make cheques or money orders payable to: Chief, Photoduplication 
Service, Library of Congress. 
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Figure 6—Mass transfer coefficient in agitated vessel vs. 
Reynolds number at 25°C. 


transfer data for engineering purposes. The general form of 


this correlation is 
SR atSeW CRO ec 6. och ee See ae (11) 


where a, / and ¢ are constants to be determined from the data 
and the definition of the dimensionless groups is as follows: 
kyl Me rd? L 

Sha =} Sow ey Be 

D prD Me 


where / = vessel diameter; d = rotor diameter; r = revolutions 
per second. 

The present data were correlated with 4 taken as 0.25, 
0.30, 0.333, 0.356, 0.42 and 0.50 respectively. The data up 
to Re = 60,000 only were considered. 

The exponents 0.25, 0.42 and 0.50 showed a wide scattering 
of the results which indicated that these exponents were unsuit- 
able for the correlation. A much better representation was 
obtained with the values 0.30, 0.333 and 0.356. The equation 
of the best straight line was obtained by the method of least 
squares in these three cases. The mean deviation of Sh number 
for the data is given below: 


Equation Mean Deviation of Sh 


(Sh) = 0.651 (Sc)®3 (Re)? -8 + 6.1% ~ sR 
(Sh) = 0.402 (Sc)®-333 (Re)0-6 + 5.8%; - (13) 
(Sh) = 0.309 (Sc)9-356 (Re)0-66 + 6.7% (14) 
Thus, the exponent b = 0.333 was thought to be the most 


suitable one to correlate the present data with the minimum 
scattering. 

\s shown by Potter®® and Friedlander and Litt” in the 
case of laminar boundary layer flows on a flat plate, theoretical 
considerations predict the Sherwood number for mass transfer 
to be given in the form 


(Sh) = 0.332 (Sc) (Re)q!!? (15) 


\ comparison of the — Equation (13) with the theoreti- 
cally derived Equation (15) of the laminar boundary layer 
theory, suggests that the two disagreed on the effect of the 
Re number on mass transfer rates, which in turn indicates that 
the /aminar boundary layer concept may not be true in the case 
of a fully baffled agitated vessel; this probably should be expected. 
\n inquiry could be then made into the predictions of the turbulent 
boundary layer theory as discussed elsewhere in terms of Equa- 
tion (7), namely, 


Sc) = Wf/2... ee) 


Based on a chart presented by Sherwood and Ryan“ for - 
evaluation of @(y*, Sc) as a function of Sc up to y* = 200, 1 
may be shown that in the range 1000 < Se < 10,000 
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Figure 7 — Reaction factor data compared with theories 
benzoic acid-NaOH at 25°C. 


(200, Sc) ~ SO". ...2.5. (16) 
Hence, Equation (7) may be written in the form 
ki(ScP?8~ — Vf72. Voge (17) 
Multiplying through by //D and rearranging 
kl l ats 
—m~ { — } (Sc)"3 (Re)a Vf/2 (18) 
D d _— 
where / = vessel diameter, and d = stirrer diameter. If //d is 


a constant, then Equation ('* could be written in a form similar 
to Equation (13), namely, 


(Sh) = a(Sc)!3 (Re)a VF/2 (19) 
where a is a numerical proportionality constant. 


Again, regarding the Sc number effect the agreement between 
theory and experiment is striking. Unfortunately, the effects 
of the Reynolds number cannot be compared because there are 
no skin friction factor data available for a fully-baflled agitated 
vessel. However, from what is known“ about fin fluid flow 
under turbulent conditions, a Re number effect such as demon- 
strated by the exponent 0.65 Equation (13) may not be 
unlikely. 

(b) Mass transfer with chemical reaction in the agitated 
vessel. 

Rates of neutralization measured were tabulated in ‘Table 2. 
Mass transfer coefficients, k,, were calculated by taking into 
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Figure 9—Comparison of reaction factor data with theories 
benzoic acid-KOH, 25°C. 
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Figure 8—Reaction factors measured by Friedlander and 
Litt compared with theories: benzoic acid-NaOH, 25°C, 


account changes in solubility of the acids due to the presence 
of the reaction product at the interface (see Figure 1). The 
corresponding k,° values were obtained from Figure 6 and 
were corrected for any change in the physical properties of 
the systems by Equation (13). Diffusion coefficients of all 
reactants used were measured in solutions of the reaction 
products and reported elsewhere“. Reaction factors, kz/kz°, 
were tabulated in Table 2 and plotted in Figures 7, 9 and 11 
for the three systems investigated here. 


1. System benzoic acid-sodium hydroxide. 

Results were plotted in Figure 7. Reproducibility of duplicate 
runs was on the average + 6.6% and on the maximum + 9%, 
he reaction factor tended to unity at low C4; values as 
would be expected. At high caustic concentrations (Cz >0.13N) 
the dissolution rate was so fast that the surface area could not 
be maintained constant and wide scattering occurred. Data 
reported by Sherwood and Ryan“ for dissolution of benzoic 
acid cylinders in sodium hy droxide solutions were plotted in 
Figure 7, too, for comparison with the present data. Very good 
agreement was observed between the two sets of data, not- 
withstanding the difference in the apparatus and techniques used. 

The data were compared with the predictions of the theoreti- 
cal models discussed previously. It may be noted that the 
stagnant-film theory line (Hatta) was calculated by Equation 
(3), with the diffusivities Dy and D, taken without ionic-effect 
considerations. This was taken care of in the top curve of 
Figure 7 which was calculated according to the ion-diffusion 
theory of Sherwood and Wei, 


The bottom curve of Figure 7 


‘ was calculated by Equation 
(10) of the boundary layer theory, 


namely, 
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Figure 10—Reaction factors measured by Friedlander and 
Litt compared with theories benzoic acid-KOH, 25°C. 
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Figure 11—Comparison with theories of reaction factors 


for salicylic acid-NaOH, 25°C. 
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The solid curve in Figure 7 was calculated by Equation (10), 
too, but allowance was made for the ionic effect on Ds and the 
viscosity effect on D, from a calculation of the sodium benzoate 
concentrations at the interface and parallel diffusivity data. 


gin the region studied, the data appeared to be better repre- 
waned by the solid curve (boundary layer theory with ionic 
effects on diffusivity) rather than by any other of the theoretical 
lines. It would be legitimate, therefore, to conclude that, in 
this particular case, the data give support to the mechanism 
suggested by the boundary layer theory, and that the experi- 
mental results justified the assumptions made in the derivation 
of Equation (10). 


Litt and Friedlander®® reported an experimental study of 
and cylinders of 


reaction factors measured for flat surfaces 
benzoic acid supported in a flow tunnel. 

Their data were recalculated and plotted in Figure 8, along 
with the predictions of the theories which were slightly different 
from Fi igure 7, since the salt concentration C4, was pr actically 
zero. 

In the region 0 < Cs/Ca; < 1.5, the data were in good 
agreement with the boundary layer theory corrected for tonic 
effects. At higher concentrations, however, the data seemed 
to be better represented by the stagnant-film line, which was 
calculated with no ionic effects taken into account. The data 
were rather limited in number, and reproducibility should be 
expected to be of the order of at least + 10% in the higher 
concentration range. Litt and Friedlander operated their 
apparatus expecting to be in the laminar boundary layer flow 
region (Re, < 50,000). As they concluded from their phy sic * 
dissolution experiments, however, turbulence was induced i 
their apparatus reflecting in a substantial increase of the aaa. 
tion rate. The treatment of the reaction factor data by Litt and 
Friedlander®?® was rather oversimplified. “hey assumed that 
Se, and Seg remain constant over the whole range of Cs/C4; 
studied. 

Actually, both Sc, and Scy vary considerably with C3/C4; 
if one considers ionic and viscosity effects. By the viscosity 
effect alone (because of the —— of benzoate at the interface), 
Se, attains a value of 1460 at C,/C4; 5, hence Seg is calculated 
to be 490 instead of 296 pide by L itt and Friedlander®®. 
For a further discussion of these data the reader is referred to 
Reference (24). 


2. System benzoic acid-potassium hydroxide. 


Data of this work listed in Table 2 were plotted in Figure 9 


and compared with the theories. Reproducibility of duplicate 
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Figure 12—Reaction factors for diaphragm cell I compared 
with calculated values. 


runs was on the average + 6.8% and on the maximum + 13.5% 
at high Cy concentrations. The data were below the Hatta line 
in contrast to the data of Figure 7 which were above. The 
deviation of the data from the boundary layer theory line with 
ionic correction was, on the average, + 3.8% with a maximum 
of + 11%. Similar, but limited, data reported by Litt and 
Friedlander®°® were plotted in Figure 10. These data were 
considerably below the Hatta line, too, and they seemed to be 
about half way between the two boundary-layer theory lines. 


3. System salicylic acid-sodium hydroxide. 

Data of this work were listed in Table 2 and aie in 
Figure It. Reproducibility was on the average + 8 . The 
lower precision was due to the smaller solubility fe salicy lic 
acid and hence to greater relative analytical errors. The data 
shown in Figure 11 do not agree with any of the theoretical 
models. The boundary layer theory line corrected for ionic 
effects was about 15% below the data and almost parallel to 
them. There are no other data reported in the literature for 
this system for comparison. Some data reported by Hsu®” 
were considered as unreliable because of poor self-consistency 
possibly due to the unsteady-state method of measurement. 


With regard to the data shown in Figure 11 there is still 
some doubt concerning the solubility data‘ of salicylic acid in 
salicylate solutions whence C,; was obtained. These 
were calculated from activity coefhcients 
»btained from sane reported at 20. rc. and the pea 
value in water taken as 0.0159M., at 25°C. Solubility measure 
ments at 25°C. made in this work from super- -saturated sohihans 
of salicylic acid in sodium — late solutions agreed with the 
predicted values within 2 ~ 3%, but measurements made from 
under-saturated solutions were considerably lower. 


sodium 
solubilities at 25°C. 


(c) Reaction-factor measurements in diaphragm cells. 


Certainly, in an agitated vessel with baffles the hydrodynam- 
ics 1s complicated and at the present not well understood. 
However, in a well stirred, diaphragm-type diffusion cell, there 
is unanimous agreement“? 5:2*:25.2 that the process of diffusion 
takes place ina truly stagnant layer of liquid with no complicating 
convention effects. 

In this work, reaction factors were measured in a diaphragm 
cell in an effort to test the validity of the stagnant-film, ion- 
diffusion theory. Experimental results obtained for the system 
benzoic acid-sodium hydroxide in two cells were listed in 
able 3. The reaction factors were determined using Equation 

20) the derivation of which is given in Reference (24), 


ky J3(1 Vi + 1/V2) 


= 20) 
Ry 2.303 DaCaK ‘ 
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Calculated values of reaction factors were obtained from 
Equation (21) below 


ky Dz Cp 
—— = 1 + Saute ewan eS aan 
ky Dy Cy 
where C, is the arithmetic mean concentration of sodium 
hydroxide at the beginning and the end of the experiment. 

The data were plotted in Figure 12, and compared with the 
calculated values. The concentrations of the reaction product 
in the diaphragm were estimated” and values of Dg and Dy 
were obtained from parallel diffusivity experiments. The 
calculated values were within + 7.4% of the experimental 
results, while the Hatta line was about 11% below the data. 
It was concluded that there is favourable agreement with the 
stagnant-film, ion-diffusion theory of mass transfer. 

Another important remark to be made is that, in these 
experiments, the “collision effect” between the counter-diffusing 
columns of OH™ and benzoate ions, if any, ought to be very 
small in this case. Perhaps, the effect of this phenomenon 
might be more evident at higher reaction product concentrations. 


Conclusions 

On the basis of the experimental data obtained in this investi- 
gation for mass transfer coefficients without chemical reaction 
for solid-liquid systems in baffled agitated vessels, it was 
established that the Sherwood number depends on the 4 
power of the Schmidt number and the 0.65 power of the stirrer 
Reynolds number. The present data are about 10% lower than 
similar data reported by Johnson and Huang“*. 


The exponent 4 on the Se group is significant in that it is 
in good agreement with the prediction of the boundary layer 
theory of mass transfer. The exponent 0.65 on the Re group 
is in disagreement with the /aminar boundary layer theory but 
it was thought to be in qualitative agreement with the turbulent 
boundary layer concept. 


Reaction factors in agitated vessels measured for solid-liquid 
neutralization systems were found to depend strongly on the 
concentration ratio of the reactants, C,/C4;, in a manner quanti- 
tatively predicted by the boundary layer theory corrected for 
ionic effects on the diffusivities of the reacting species. The 
data were in agreement with similar data from rotating dissolving 
cylinders, a fact which indicated that the geometry of the 
apparatus is insignificant. 

Reaction factors in diaphragm-type diffusion cells were 
obtained and were found to quantitatively support the theory 
of ion diffusion in a stagnant film. 
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Nomenclature 

a,b,c = Constants 

5 = Concentration, g. moles cm.~‘ or g. equiv. cm 

D = Diffusion coefficient, cm.? sec. 

d = Stirrer diameter, cm. 

Ep = Eddy diffusivity, cm.* sec.~! 

f = Skin friction factor. 

a = Solute mean flux, g. moles. cm.* sec.-! or g. equiv. 
cm,~? sec,~! 

ki = Mass transfer coefficient with chemical reaction, cm. 
sec,~! 

ky? = Mass transfer coefficient without chemical reaction, 
em, sec.~! 

K = Diffusion cell constant, cem.~? 

l = Agitated vessel diameter, cm. 

N = Rate of mass transfer, g. moles. cm.~? sec.~! or g. 
equiv. cm.~* sec.~! 


r Revolutions of stirrer per second 
(Re)j¢ = Stirrer Reynolds number, rd*p/p 
Rate of surface renewal, sec.~! 


ll 


ll 
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Sc = Schmidt number, u/pD 

Sh = Sherwood number, kz°//D 

St = Stanton number, k7°/i 

T = Temperature, °K. 

ii = Mean fluid velocity, cm. sec.~ 

Vi,V2 = Volumes of compartments 1 and 2, cm.’ 

x = Distance along diffusion path, cm. 

XL = Thickness of diffusion layer, cm. 

y = Distance from wall, cm. 

y* = Dimensionless distance from wall defined by Eq. (9) 


Greek Letters 


a = Constant 

€ = Eddy kinematic viscosity, cm.? sec.~! 
0 = Time, sec. 

be = Fluid viscosity, poises or g. cm. a. —" 
v = Kinematic viscosity, H/p, cm? ‘ 
p = Fluid density, g. cm.~ 

Te = Shear stress at the wall, dyn. cm.~? 
d( ) = Denotes function 

Subscripts 

A,B = Solutes A,B 

D = Diffusion 

l = Interface 

L = Liquid 

w = Wall 
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High-Speed Photography in the Study 
of Two-Phase Flow’ 


R. A. S 


Some of the possible applications of high-speed 
motion picture photography in the study of two-phase, 
two-fluid flows are outlined. These include measure- 
ments of ripple or wave velocities in annular flow 
and bubble velocities in bubble and slug flows. 
Qualitatively, the motion pictures permit better 
understanding of the interactions between the two 
phases. 


The equipment being used in the present photo- 
graphic study of vertical flow of oil-water mixtures is 
described. A 1l6-mm. “Fastax” high-speed motion 
picture camera having speeds to 8000 frames/sec. is 
being used in these studies. Auxiliary equipment 
includes a high intensity light source and a timing 
pulse generator. 


In the oil-water system being studied, bubble 
groups and bubble clusters are formed which rise 
at approximately constant velocities. Very little effect 
of bubble size on the velocity of the bubble was 
noted except in a few instances where the bubbles 
appeared to be close to the wall. Preliminary results 
of the effects of superficial oil and water velocities 
on the bubble velocity are presented. 


Some other photographic techniques which might 
find application in the analysis of two-phase flows 
are discussed. 


[' has been apparent for some time that a need exists for a 


more fundamental approach to the analysis of the flow of 


two-phase mixtures. While numerous empirical correlations 
have been proposed within the last few years“:® no great 
improv ement in the understanding or the accuracy of prediction 
of flow characteristics has resulted. There have, however, 
been a few notable attempts to approach the analysis of individual 
flow patterns from a fundamental basis“. Mechanistic 
analyses such as these are promising and should eventually be 
more productive than the strictly empirical approach impor- 
tant though that a is. 


\t the University of Alberta, high- speed motion pictures 
are being taken to serve as a guide in the analysis of the flow 
of two- phase systems. Motion picture photography is useful 
for this purpose since the interactions between phases can be 
observed directly and accurate measurements of the flow geom- 


‘Manuscript received January 17; accepted May 5, 1961. 

*Ph.D. candidate, Department of Chemical and Petroleum Engineering, 
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‘Professor of Chemical Engineering, Dean of Faculty of Engineering, 
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etry and the phase velocities may be obtained. To date, most 
of the tests have been concerned with the vertical flow of oil- 
water mixtures in the bubble and slug flow regions. 


In bubble flow, high-speed motion pictures can yield i 
formation on individual bubble and bubble cluster caedaian 
and their paths of rise, distribution of bubble diameters, in situ 
volume ratio of the phases, interactions between bubbles and 
interactions between the phases. The latter would include 
oscillation of the bubbles if the frequency was sufficiently low 
to be observable. 

High speed motion picture photography should be equally 
productive of useful information in the other flow regions. A 
survey of the measurable characteristics would include bubble 
velocity and shape in slug flow and ripple or wave velocity, 
amplitude and period in annular flow, as well as in situ volume 
ratios of the phases. In other regions, such as the froth flow 
region, so little is known about the flow that possible measure- 
ments can only be surmised until some tests are made. For 
some of these cases, also, special methods of lighting may be 
required to outline the flow structure. 

The above measurements are some of those obtainable from 
the motion pictures. In all of these, the main purpose is to 
obtain the variation of the observed characteristics with the 
flow rates and the properties of the two phases and to learn 
from this more of the mechanisms important in the different 
regions. 


The Bubble Flow Region—Theoretical Considerations 


Consider a mixture of two phases one continuous, the 
other dispersed as bubbles entering a vertical section of tube 
in time t*. Let the superficial velocity of the continuous (refer- 
ence) phase be Ve and that of the dispersed (second) phase be 
Vs. The flow is in the bubble or slug flow region with N 
bubbles entering the section in the time f* and having individual 
volumes 7; and instantaneous velocities V si. Assume that none 
of the bubbles introduced in the time interval overtake bubbles 


introduced before the time t = 0. Att = f*, the bubbles will 
extend over some length L where 
= fo'* Vopdt.. (1) 


and bubble p has moved furthest up the tube. The volume of 


second phase which has entered is 
N 
dr = VsAd*. s ‘ os ‘ (2) 
Q 

and hence the volume of continuous phase is 


LA, — VsAd* = Ajfol*Vigdt — VsAd®... (3) 
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Figure 1—Schematic diagram of photographic equipment. 


If ¢* is sufficiently great, the volume ratio of second phase to 
reference phase will be characteristic of the flow and the holdup 
ratio (1) may be written 


J 's J ’R 


Hr = —_———* 
- (At of* Vopdt = 


A,Vst* 


A, Vst* ) 


(* Vi pdt _ Vst* 
=— = Soutien , (4) 


Vrt* 
For the special case where all bubbles are moving with the same 
velocity, V;, Equation (4) becomes 


ge a 


The average density of the mixture can be written as 


Vst* Sut Vepdt — Vst* 


= ps a ; (¢ 
wie Vopdt - fol* Vi pdt " 


Po 
or, for constant bubble velocity 


Vs i acl i 
Ps = Ps y, + PR — 


=p: - (pr = Os) 











Figure 2—Photographic equipment and test section. 
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‘The pressure drop over a section of vertical tube is composed 
of a hydrostatic head component and a frictional component. 
For the bubble and slug flow regions, the hydrostatic compo- 
nent of the unit pressure drop is equal to the av crage densi- 
ty of the mixture in the tube and 


AP 


hee oe ome 
rw (8) 


= Pa. 


In the bubble flow region, the frictional pressure drop should be 
very closely approximated by the pressure drop calculated for 
reference fluid flowing alone at the mixture velocity, Vm. 


AP 2fprVn 
=f = a pens S 6: o S68 . eee eo ewe 9 
ie. — gD (9) 
where Vm = Vr + Vs 


and f depends upon a Reynolds number calculated as 


D\ ‘mPR 
MR 


Re = 


The calculation for slug flow is the same in principle although 
some differences exist due to the restricted cross-section for 
flow of the reference phase. 


Experimental Equipment 


The experimental equipment is shown schematically in 
Figure | and photographs of the working level of the eqpmnen 
— and of the control equipment are given in Figures 2 and 3. 

\ description of the flow system being used for the present 
study of oil-water mixtures has been given previously“. The 
test section is a cellulose acetate butyrate tube having an inside 
diameter of 1.038-in. The flow tube is 37-ft. long and the 


Figure 3—Control equipment. 
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Figure 4—Viewing section. 


viewing section 1s 23-ft. above the mixing tee. The test section 
for pressure drop is 28.37-ft. long and the lower tap is about 
8-ft. above the mixing tee. Properties of the oil and water 
are given in Table 1. 

The tube has been inscribed in the viewing section with 
lines 14, 1, 3 and 6-in. from a reference line. This section of 
tube is enclosed in a view-box shown in Figure 4 consisting of 
a two-foot length of 3-in. 1.D. Lucite tube which has been cut 
lengthwise along a chord and faced with '¢-in. Lucite sheet. 
End pieces are Lucite discs which were machined to fit over the 
flow tube and centre it in the view-box. A liquid seal around 
the flow reget is obtained with O-rings set in grooves in the end 
pieces. A drain hole with tap at the bottom of the view box 
and air holes at the top permit addition and removal of liquid 
around the flow tube. 


Photographic equipment was purchased as a unit prepared 
by the Wollensak Opt ical Company and consists of camera, 
control or “Goose’’ unit, high intensity Xenon lamp unit and 


100/1000 cycles per second timing pulse generator. The 
oa = 
we - is 





Figure 5—Projection equipment used in analysis. 
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TABLE 1 
Puy SICAL . PROPE RTIES OF THE LIQUIDS 









Water 
Densitv 
Viscosit y 


62.3 Ib./ft.8 at 72.0°F. 
0.936 cp. at 72.0°F. 





Oil (Marcol GX) 
Density 
Viscosity 


53.00 Ib. /ft.? at 72.0°F. 
20.09 cp. at 72.0°F. 


Interfacial Tension 50.34 dynes/cm. at 77.0°F. 





camera is a Fastax WF 3 model which uses 100-ft. reels of 
16 mm. film. In this camera, which has a maximum repetition 
rate of about 8000 frames/sec., the image is stabilized on the 
continuously moving film by a prism w hich rotates in synchro- 
nism with the film motion. Different film speeds are obtained 
by adjusting the voltage applied to the universal motors of the 
camera and the repetition rate at any point is indicated by timing 
pulses registered on the film during operation. A prism arrange- 
ment permits through-the-lens focussing. The high intensity 
light source is a discharge lamp which can be pulsed at three 
power levels or used as a semi-continuous source at a fourth, 
lower, power level. Direct current for the lamp is supplied by 
a bank of six 6-volt automotive batteries having a fairly high 
current rating. 

The pulse lamp was positioned directly behind the viewing 
section as shown in Figures 2 and 4 since preliminary tests 
indicated that better contrast could be obtained with back 
lighting. It was necessary to shield the camera from direct 
light from the lamp (which might cause fog ogging of the film) by 
installing light shields on each side of the view box. 

Kodak Tri-X negative film (Type 7233) has been used for 
most of the tests presented here although some tests have been 
made with DuPont Superior 4 Negative Film (Type 928A). 
lhe processing is quite simple for either of these films. 

The film was analyzed by projecting it on a small rear- 
projection screen as shown in Figure 5. The projector used 
was a Wollensak Motion Analysis machine which permits 
single frame projection, has fairly good registration and is 
equipped with a frame counter. The measurements of the 
flow geometry and velocities were aided by superimposing a 
grid on the projected image. This grid was a photographic 
copy on Mylar film of a “sheet of standard arithmetic graph 
paper (10 X 10 to the !4-in.) and was placed directly against 
the rear projection screen. 


Experimental Procedure 


For a given position of the camera with respect to the 
viewing section, the combination of lens and extension tube 
used was determined by the field of view desired which, 1 
this case, was the length of viewing section to be observed. 
This was dependent, in turn, on the information desired from 
the film. If the viewing length is long, an individual flow 
configuration can be observed to determine the chan; ges in its 
behavior with time (or position) and its average “velocity. 
Shorter viewing lengths (greater magnification) are better 
suited to the determination of geometrical characteristics of 
individual configurations and instantaneous velocities. 

Ihe next factor to be considered was the picture repetition 
rate desired. A lower limit on the repetition rate is set by the 
velocity of the pattern being observed and the resolution of the 
optical system being used. Stated simply, if the subject moves 
a distance greater than the resolution limit of the optical system 
during the period of exposure the image will appear blurred. 
(A good discussion on the requirements for “‘stopping’’ motion 
by high speed motion picture photography is given by Hyzer™.) 
Phe choice of film speed is also determined, therefore, by the 
purpose of the test. If the study is directed towards a study of 
bubble vibration or oscillation rather than towards the deter- 
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TABLE 


2 


_ EXPERIMENTAL ConpDITIONS 


Oil Rate 


























Test No. | Ma ae — Pressure Drop Holdup Ratio 
| ft. H.O/ft. Hr 
ft.3/sec. ft./sec.@ ft.3/sec. ft./sec. 
18 0.00586 1.0 0.00127 0.2165 0.990) 1.45) 
19 0.00586 1.0 0.000434 0.0739 1.002 Las 
26 0.00586 1.0 0.000339 0.0577 | 1.00 1.67 
27 0.00188 0.32 0.000434 0.0739 0.990 2.20 
28 0.0106 1.8 0.000434 0.0739 1.014 1.352 
29 0.00586 1.0 0.000902 0.154 0.998) 1.50 
31 0.00586 1.0 0.00507 0.865 | 0.960 1.44 
34 0.00586 1.0 0.00264 0.449 | 0.975) 1.42) 
35 0.00586 1.0 0.00507 | 0.865 0.960 1.44 
@Supe rfic ial lineal velocity 
) Interpolated from Reference (9) 
TABLE 3 
PHOTOGRAPHIC CONDITIONS 
a ; Distance eh Nominal Lens Aperature 
— | Camera to Tube Lens Lighting Repetition Rate f 
18 8 ft. 2 in. 153.6 mm. 2.5 kw pulse 1500 32 
19 8 ft. 2 in. 153.6 mm, 2.5 kw pulse 1500 32 
26 29 3% in. 50 mm. Simmer 1500 il 
27 | 29 3% in. 50 mm. | Simmer 1500 11 
28 29 3% in. 50 mm. Simmer 1500 a6 
29 29 % in. 50 mm, Simmer 1500 5.6 
31 8 ft. 2 in. 153.6 mm. Simmer 1500 6.5 
34 8 ft. 2 in. 153.6 mm. Simmer 1500 6 
35 8 ft. 2 in. 153.6 mm. Simmer 1500 6.5 
| | 


(@) Pulse lamp 1 foot behind and slightly above viewing section. 


mination of average bubble velocities, higher repetition rates 
will be required because of the greater velocity of the observed 
phenomenon. 

Knowing the repetition rate, the lens setting for a given 
lighting arrangement may be determined. Some ‘difficulty was 
experienced with the choice of Jens settings for the back lighting 1g 
arrangements used because the incident light meter gave settings 
which were too low. Experience should indicate, however, the 
correction which must be applied to the meter reading. 

Experimental flow conditions for the tests reported here are 
given in Table 2, while the photographic conditions are noted 
in Table 3. As the main interest in these tests was in average 
velocities of bubbles and bubble groups, a relatively low repeti- 
tion rate was used. 


Experimental Results 


The lighting method used gave fairly good contrast and the 
bubble motion is easily observable. The bubble image is rather 
hazy and the reason for this has not yet been determined. This 
lack of definition does not appear serious except in the size 
determination of small bubbles where the percentage error may 
be large. Velocity measurements are quite simply made and 
there appears to be very little possibility of judgment errors in 
the determination. 

These initial investigations with high-speed motion picture 
techniques have been made in the bubble and slug flow regions 
for a system consisting of a white oil and water. The set of 
tests comprising the present investigation includes six tests at 
the same oil velocity and varying water velocity and an additional 
two tests at one of the water velocities and different oil velocities. 
One check run was also obtained at a high oil flow rate. The 
water and oil flow rates are given in Table 2 along with the 
pressure drop and hold-up values for the test conditions as 
obtained from earlier results“ 

The most notable flow characteristic obtained from the films 
was the relative constancy of the velocity of rise of bubble 
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groups. The size of the individual bubble has no observable 
effect on its velocity of rise when in a bubble group. This 
effect, which might be expected when the bubbles are clustered 
together would not necessarily be expected for more open 
groupings of bubbles where the individual bubbles might be 
expected to rise in relation to their sizes. The result indicates 
that it may be more appropriate in some cases to treat bubble 
flow as a type of slug flow rather than as a modified single 
bubble flow. 


Figure 6 shows the change in bubble size and grouping for 
superficial water velocities of 0.32, 1.0 and 1.8-ft./sec. with a 
superficial oil velocity of 0.0739-ft./sec. Similarly, the flow 
patterns for six superficial velocities at a superficial water 
velocity of 1.0-ft./sec. are shown in Figure In each case 
the leading end of a group of bubbles is shown. Figure 7 
indicates the development of slug flow as the flow rate of the 
second phase (oil) 1s increased at a constant flow rate of reference 
phase (water). The highest oil flow rate is well within the slug 
How region and approaching froth flow. The transition between 
regimes | and Il occurs at a superficial oil velocity of about 
1.25-ft./sec. 

The bubble velocities obtained in the present investigation 
are given in Column 2 of Table 4. Columns 3 and 4 show the 
comparison between the holdup ratio obtained directly“) and 
that calculated from the bubble velocity with Equation(5). The 
average density of the mixture in the tube, obtained from 
quation (7), is given in column 5 while columns 6 and 7 show 
the comparison between the experimental pressure drop and 
that calculated from Equation (9). 

Figure 8 indicates the effect of superficial water velocity on 
the bubble velocity for a superficial oil velocity of 0.07 39-ft./sec. 
lhe dashed line indicates the mixture velocity. The slip velocity, 
defined as the difference between the bubble velocity and the 


mixture velocity, increases from 0.39-ft./sec. at a mixture 
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Figure 6—Flow patterns for a superficial oil velocity, V,, 
of 0.0739-ft./sec. 
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velocity of 0.39-ft./sec. to 0.53-ft./sec. at a mixture velocity of 
1.87-ft./sec. 

The effect of the superficial oil velocity on the bubble 
velocity is shown in Figure 9. In this figure also the mixture 
velocity has been given and, in addition, the variation of slip 
velocity with superficial oil velocity is indicated. In the region 
of bubble flow, which extends to an oil velocity of about 
0.40-ft./sec., the bubble velocity is relatively insensitive to the 
oil rate. In the slug flow region, however, the influence of the 
oil rate on the bubble velocity becomes quite significant. The 
lower curve shown in Figure 9 indicates the variation in slip 
velocity with superficial oil velocity. It is interesting to note 
that this relationship passes through a minimum at about the 
bubble-slug transition. 

Figure 10 shows a first estimate of bubble size distribution 
for a “superficial oil velocity of 0.154-ft./sec. at a superficial 
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Figure 7—Flow patterns for a superficial water velocity, V ;, 
of 1.0-ft./sec. 
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Figure 9—Effect of oil velocity on bubble velocity. 


basis of 
the fre- 


both on the basis of the number of bubbles and on the 
the volume of bubbles. While the size distribution (on 
quency basis) should be substantially correct for the smaller 
diameters, there appears to be considerable error in the estima- 
tion of diameters of large bubbles. The flow rate of oil obtained 
from the size distribution given is approximately 30% low. 
This large discrepancy is not too surprising in view of the 
difficulty involved in estimating from a two dimensional view 


water velocity of 1.0-ft./sec. The size distribution is shown the equivalent diameters of the highly distorted, large bubbles, 
TABLE 4 
EXPERIMENTAL RESULTS 
Holdup Ratio 
Bubble Average Pressure Calculated 
Test No. Velocity Density” Drop Pressure Drop 
ft./sec. Observed From Bubble | Ib. /ft.3 Ib. /ft.2/ft. Ib. /ft.2/ft. 
Directly Velocity™ | 
18 1.60 1.45 1.38 61.1 61.7 61.6 
19 1.55 1.73 1.48 61.8 62.5 62.2 
26 55 1.67 1.61 61.9 62.3 62.3 
27 0.780 | 2.20 2.20 61.3 61.7 61.4 
28 2.40 1.52 | 1.29 62.0 63.2 63.1 
29 ‘37 1.50 1.42 61.4 62.2 61.9 
31 2.40 1.44 1.55 59.0 59.8 60.1 
34 | 1.76 1.42 1.31 59.9 60.75 60.6 
35 2.30 | 1.44 1.44 58.8 59.8 59.9 


“From Equation (5) 
From Equation (7) 
©From Equation (9) 
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Discussion 


High- speed motion picture photography appears to have 
certain advantages in the study of two-phase flows. The most 
prominent of these is the possibility of obtaining simultaneously 
the geometry and behaviour of individual flow units. The 
analysis of the motion pictures may be carried out quite readily 
if the projection equipment used has a frame counter and can 
operate in single frame operation. The use of a rear projection 
screen with a superimposed grid appears to have certain advan- 
tages for the analysis of film if sufficient accuracy can be obtained. 


In earlier investigations” it has been assumed that the 
entrance section had no effect on the resulting flow which was 
a function only of the flow rates of the phases, the phase 
properties and the tube size. This assumption was based upon 
preliminary tests made with various entrance or mixing sections 
with the air-water system. It now appears possible that the 
entrance section may be important in some cases, particularly 
for immiscible liquid systems at low flow rates. It is difficult 
to imagine the grouping and clustering of bubbles observed at 
low flow rates arising if the oil were injected axially through a 
smaller tube. At higher flow rates, however, when the rates of 
shear in the mixing tee are high, the flow pattern undoubtedly 
is independent of the entrance section. 

The good comparison of observed and calculated holdup 
ratios and pressure drops indicates that the assumptions made 
in the calculations are reasonable (and the bubble flow can, in 
this case, be treated as a modified slug flow). Since the frictional 
component is such a small fraction of the total pressure drop, the 
analysis of the frictional losses could be considerably in error. 
The theory presented should apply to any periodic flow if a 
characteristic velocity exists for the flow units. It is interesting 
to note that in all but two cases the observed holdup ratio is 
greater than that calculated from the bubble velocity. As the 
observed holdup ratio is the average for the entire vertical flow 
section, including the entrance section, this behavior is in agree- 
ment with the results of Moissis and Griffith® who predict 
significantly higher hol dup ratios in the entrance section for 
slug flow. No quantitative comparison is possible, however, 
since their data apply to gas-liquid flow. 


‘There are a number of possible modifications of the equip- 
ment which might make the data more valuable. For some 
investigations a measure of pressure or pressure gradient ob- 
tained along with the motion picture record of the flow might 
prove valuable. This could be obtained most easily by installing 
pressure transducers in the test section being photographed and 
connecting the output leads to an oscilloscope mounted beside 
the test section. There are also cameras available which will 
record an oscilloscope trace on the film through a separate 


optical system. 
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(7) Hyzer, W. 


Another modification which is being considered at present 


is the installation of a mirror system to give a a side view of the 
test section. This would be a great help in estimating bubble 
volumes since, as mentioned above, the bubbles are frequently 
not symmetrical. 

As has been mentioned above, special methods of lighting 
may be required to outline the flow structure in some cases. 
This referred particularly to the use of a transverse plane of 
light to outline the surface of the water film on the tube wall 
for slug and annular flows. It may also be possible to use this 
technique in outlining the structure of froth flow. 


Nomenclature 
A, = Cross-sectional area of tube, ft.? 
D = Diameter of tube, ft. 
: snes (AP/ Ax) ) BD 
f = Friction factor, f = —— -—, dimensionless 
P ; 2 V2 
: ‘ ft: Ibm. 
g- = Dimensional constant, —— 
by. sec. 2 
Hr = Holdup ratio, dimensionless 
L = Length along tube, ft. 
AP , 
- Pressure gradient, lb./ft.2/ft. 
Ax 
Re = Reynolds number, dimensionless 
t = Time, sec. 
v = Bubble volume, cu. ft. 
V = Velocity, ft./sec. 
bh = Viscosity, lbm./ft. sec. 
p = Density, lb»./cu. ft. 
Subscripts 
a = Average 
b = Bubble 
p, i = Bubble designation 
= Reference (water) phase 
S = Second (oil) phase 
Superscripts 
* = Time interval 
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Determination of Ternary Liquid-Liguid 
Systems by Use of the Lever Rule’ 


E. L. HERIC? and KENNETH R. WILLIAMS: 


A method is described for the determination of 
ternary liquid-liquid systems. The approach is based 
upon a combination of the lever rule and the 
measurement of a physical property which is a 
function of composition. The method has been tested 
upon the system toluene-water-acetic acid at 25 + 
0.05° and found to yield results in agreement with 
literature values. A possible modification of the 
method—which may be necessary in systems where 
tie-lines are essentially parallel when expressed on 
the usual weight per cent basis—is described. 


a" traditional method of study of ternary liquid- liquid 
systems is that of “cloud point” titration”. For precise 
results this technique is tedious, and, moreover, there are at 
least two inherent disadvantages to the method: (a) one cannot 
explicitly fix a point on the binodal curve, but only reckon that 
it must lie between two limits, one repre senting the appearance 
and the other the disappearance of a second phase; (b) the 
appearance of the second phase may be uncertain in that its 
form may vary even in a given sy stem. The subjectivity of 
this method may be avoided by use of the procedure described 
by Prutton et al), although the ‘latter may be somewhat laborious 
in that it requires a relatively large number of known synthesized 
mixtures. An alternative to the above methods is presented 
here, which may have certain advantages; the number of weigh- 
ing operations is not large, yet it is possible to fix with the same 
series of known synthesized mixtures both the binodal curve 
and the tie-lines. The method involves use of the lever rule“ 
in combination with measurements of a suitable physical property 
which varies continuously over the binodal curve. 

‘To facilitate explanation of the present method it is described 
in terms of the system toluene-water-acetic acid; this was studied 
as a means of determining the practicality of the method. No 
tabulation of results obtained by the present method is presented, 
as the agreement only confirms the previously reported values‘. 
lo summarize, results of both studies agree to the extent that 
the average departure of the binodal points obtained by the 
present method is + 0.1 weight % from the smoothed curve 
drawn through the haa values. Present and literature 
tie-line data also are in quite good agreement. 
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The Method 


For best application of the method it is advisable to have 
an idea of the general nature of the system to be studied. This 
need consist only of the approximate location, perhaps to the 
nearest whole % by weight, of three points on the binodal 
curve, namely, the plait point and the mutual solubilities of 
the two non-consolute liquids, respectively C, De and Dz in 
Figure 1*. The cloud point method may be used to adv antage 
for this purpose; for such approximate work it is rapid and the 
objections noted above are not pertinent. The plait point was 
located by alternately adding each of the three components to 
an initial known toluene-water mixture so as to maintain about 
equal volumes of the two phases until the system abruptly 
became homogeneous. The approach to the plait point is 
accompanied by a convergence of the densities of the two 
phases; this effect serves as a useful indicator of the point’s 
proximity. 

After fixing the locations of these points, an estimate of the 
binodal curve was made, the dashed line in Figure 1. On this 
basis two series of heterogeneous mixtures were prepared. The 
first series, ‘reference mixtures”, indicated by points Rixy, 
Royy and R3yy was prepared by adding variable known amounts 
of the homogeneous mixture (X) of toluene and acetic acid to 
known amounts of heterogeneous mixtures (1’) of toluene and 
water. X and Y' were selected so that the straight line, XY, 
joining these two points the line upon w hich all mixtures 
formed from them must lie would pass through the binodal 
curve near the plait point and on that side of it which contains 
phases richer in the same non-consolute making up the major 
part of Y. The purpose of this selection, as will be evident 
below, was to allow the determination of the system over all 
ranges of concentration, including that near the plait point. 
X and ¥' should also be selected so that even without using 
overly large amounts of the components there will be sufficient 
quantities of both phases of each reference mixture for the 
procedures to be described. 


The second series of known mixtures, “scanning mixtures,’ 
indicated by points Sy, Ss and Ss, was to lie in the region toward 
the portion of the binodal curve opposite the reference mixtures. 
Scanning mixtures need not fall on a common line, but again 
there should be sufficient amounts of both phases of each 
mixture as noted above. 


\fter allowing the various thermostated mixtures to reach 
equilibrium after agitation, a sample of each phase was removed 


*Throughout this work, letter subscripts refer to various specific points. 
In this instance, U and L refer respectively to the upper and lower phases 
on the base tie-line. Number subscripts refer to the tie-lines upon which 
points fall when such identification is necessary. Letter superscripts refer 
to a coordinate of a point in terms of the component indicated. 
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Dy Y WATER > Di 


Figure 1—Reference, scanning and extrapolated mixtures. 


and the refractive index measured with a Bausch and Lomb 
precision refractometer. It was now possible to construct the 
first two plots shown in Figure 2. These are plots of mp versus 
weight % acetic acid for reference mixtures. The curves are 
two branches, each corresponding to one side of the binodal 
curve. These plots in effect fix one point on each tie-line of the 
system, namely, that point upon the line XY in Figure 1. (In 
reality only one phase of the heterogeneous region need be 
considered, but a study of both phases provides a rapid means 
of duplicate determination and appears to be worthwhile.) The 
two plots were combined with the refractive indices of scanning 


samples in constructing the third plot in Figure 2, by means of 


which a second point on each tie-line may be located. This is 
equivalent to the establishment of the slopes of the tie-lines. 
The third plot has been formulated in terms of weight % acetic 
acid on the acid-water axis versus weight % acetic acid for the 
reference mixture which falls on the same tie-line. The only 
significance of the first of these functions is that is shows where 
the tie-lines would intersect the acid-water axis if they were 
extended beyond the binodal curve. 

Establishment of points used in constructing the third plot 
is based upon the fact that any two points on the same tie-line 
yield the same pair of phases at equilibrium. Thus for mixture S; 
in Figure | the refractive indices of the two phases corresponded 
to those of a point Ryyy interpolated from the first two plots in 
Figure 2. By passing a straight line through points S; and Ryyy 
in Figure 1 it is possible to fix point E44w on the acid-water 
(AW) axis, an “extrapolated mixture”. Points Es4w and Egaw 
are similarly fixed using the pairs S; — Rsyy and Ss — Roxy. 
The extrapolation may be performed either graphically or 
algebraically. In passing, it should be noted that, after the 
third plot of Fi igure 2 has been constructed, it in turn may be 
used to fix corresponding pairs of points such as Royy — Bus 
and Riyy — Eyaw in Figure 1, thus defining additional tie-lines 
in the system. 

The basis for location of the binodal curve points is the lever 
rule, applied in a manner similar to one previously described. 
Assume a known heterogeneous mixture, /; in Figure 3. At 
equilibrium the mixture will exist as two phases, an upper 
phase /\y and a lower phase /,,, the location of which is the 
present aim. It should be noted that mixtures such as /; need 
not be in addition to those described as reference or scanning 
mixtures above, but may be either of the latter if sufficient 
amounts of them are available for the procedure to be described. 
This convenience is an advantage of the present method. 
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Figure 2—Working plot for the determination of tie-line 
slopes. 


(1) ny, versus weight % acid in reference mixture, water- 
rich phase. 

(2) ny, versus weight % acid in reference mixture, toluene- 
rich phase. 


(3) weight % acid in reference mixture versus weight 
% acid on acid-water axis. 


Determination of mp values of the phases /jy and 1), is suffi- 
cient to fix the tie-line upon which these points lie. Let a working 
mixture of known composition, J in Figure 3, now be added to 
a known amount of a sample of phase /\y, and the resulting 
mixture F», allowed to reach equilibrium. The choice of J is 
quite arbitrary except that it should cause Fy» to lie in the 
heterogeneous region. J may be a three- or two-component 
mixture, or even one of the pure components, whichever appears 
to be the most convenient. After F2 has reached equilibrium, 
measurement of mp of each phase, Foy and F2;, establishes the 
tie-line upon which F; lies. 

Suppose that a weight g 


typ Of phase Jy has been added to a 


weight ¢, of J. Let there be defined 


AcID 








WATER > 
Figure 3—lllustration of the basis for determination of 


points on the binodal curve. All ternary points are in the 
heterogeneous region of the phase diagram. 
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where the bar indicates the length of the line joining the indicated 
pair of points, it follows that 

Geet ce Me ni 

FA —_ JA 

Thus if straight lines are drawn from point J to intersect the 
tie-lines containing /,y and F», only that single line passing 
through both of the latter points will yield the values of the A 
coordinates required by Equation (3). The desired line is 
established by a trial-and-error procedure. 

Although the method may be applied in an essentially 
graphical manner, it is pertinent to express the relationships 
analytically, also. ‘This provides an approach which is not 
subject to some of the practical limitations of graphical pro- 
cedures. The tie-line upon which /,, and thus /,y and Jz, lie 
may be defined by the linear equation 

me Mh Bs. ices cas sus (4) 


where M, and B, are respectively the slope and intercept. In 
terms of quantities indentified in Figures 2, 3 


A 1AW a, TAixy 


M —_——...... 5) 
, ee me = Fe, ( 
and 
ge cic. tx tw ~tiveses dC 
In turn, for use in Equation (5) the relationships 
Line = 100 — Iiaw..... sé Gb] 
and 
yw 


PM ise (AS — TA as is ev cee wv (8) 


Xa 

are available. As Equation (8) is used repeatedly in the study 
of a system, it is pertinent to point out that all quantities except 
I*\yy are constant once X and Y have been selected for the 
system. For the tie-line containing F*,» there are the following 
relationships, aang to Equations (4) - (8) above: 


Fa, = My F¥, + Bo. (9) 
Fauyw — FAs x Y 
M, = FW, aa FW, 5 eees . ° . . (10) 
By = FAga7 ee oe ° as . (11) 
FWo,w = 100 — Fasay...... ; (12) 
yw 
FWyxy = (X4 — Fagyy) YA Tere : (13) 


If, now, Equation (14), comparable to Equation (3) above, 


iy — Ju 
p, — Je 


is combined with Equations (3), (4), and (9), there 1s obtained 
Equation (15) explicit in /4 yw: 


Mil(G — 1) (J4* — MsJ™) — GB] + MsB, . 
I4\y = — ——.... (15) 
M2 — M, 





The last equation is in terms of known parameters of the system. 
A second coordinate of the binodal point, /" vy, may be obtained 
from Equation (4). 
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Figure 3 also contains points associated with the location 
of 11, the lower phase of mixture /;, by the addition of a working 
mixture K to obtain a point F3 in a manner quite parallel to that 
just described for the location of I,y. If, corresponding to G in 
Equation (1), there is now defined 

&n, + 8x 
CO  ——. nc cdc ee bcmeiewe (16) 


uy 


where the subscripts are now pertinent to the points /,z, F3 and 
K, it can be shown that, corresponding to Equations (3) and (15), 


Gf ase (17 
“7, Ke deta (17) 
and 
M,[((G’ — 1) (K4 — M;KW) — G’B;} M;B 
ig a )¢ 3 3} + MB, -.. 8) 





M; — M, 


Point /;, may thus be located graphically by trial-and-error 
methods with ‘Equation (17), or analytically with Equation (18). 
It should be noted that there is no fixed relationship between the 
tie-lines containing respectively points F, and F3, as each is 
dependent only upon the relative amounts of Jv, /:z, J or K used. 

Aside from difficulties inherent in the use of graphical 
procedures for precise results, the limiting aspect of the present 
method is implicit in the denominators of the functions in 
Equations (15) and (18). As M2 or M3 approaches M,, that 
is, as the tie-lines approach parallelism, the uncertainty in the 
A coordinate of a binodal curve point approaches infinity. This 
may be predicted, of course, purely from the geometry of 
parallel lines. It can be shown, however, from the relation 
between an angle and its tangent, that a difference in slope of 
only six degrees between the two tie-lines used in establishing 
a binodal point is sufficient to reduce the uncertainty in the 
latter to about + 0.03 weight %. The uncertainty decreases 
as the slopes become more unlike, becoming, for example, about 
+ 0.01 w eight % when the difference is 17 degrees. In stating 
these uncertainties some assumptions have been made concerning 
the precision of measurements. First, the compositions of 
synthesized mixtures should be known to the nearest 0.01 
weight %. ‘This requirement is quite readily satisfied even 
in the preparation of mixtures of pre- selected composition, Y 
in Figure 2, by the use of microliter pipets for final adjustments. 
Second, weights of Iv, Ii, J or K in Figure 3 should be 
determined to | part in 20,000. These requirements concerning 
weights and compositions mean that the volume of the vapor 
phase must be kept small, or that it must be possible to correct 
for the weights of the components in the vapor phase. Finally 
a G value of no more than about 1.5 should be used if extrapola- 
tion errors inherent to such a method are not to become overly 
large. 

It is evident that, even though the method as described 
would be applicable to many systems, its usefulness would be 
increased if it could be applied to systems where tie-lines show 
only a small departure from parallelism on the usual weight % 
plot. While this complication was not a factor in the system 
described here, the consequences have been considered in antici- 
pation of such a problem. A possible solution to this problem 
lies in a transformation of concentration variables from weight 
% to mole % at a certain stage in the procedure. Such a 
transformation will produce a spectrum of tie-line slopes in 
these systems when there is a significant difference between 
the molecular weights of the non-consolute components. 


In modifying the method as previously described on the 
basis of weight %, no change need be made until the aspects 
of Figure 3 are considered. Up to that point use of weight % 
is quite acceptable, and undoubtedly more convenient. Consider 
Figure 3 again, but now visualize that each designation refers 
to mole %. To distinguish between the previous interpretation 
of Figure 3 and that now to be assumed, suppose that all of the 
labelled points are altered from upper to lower case letters — 
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the latter to indicate the same quantities as before but now in 
terms of mole rather than weight %. Thus, for example, the 
equation of the tie-line containing i; may be written as 


14, = mi; +h... en. . ia 8 6 (19) 


where mm and /, are respectively the slope and ianecpest: Other 
equations similar to those stated in terms of weight % may be 
listed, and the present modification of binodal point determination 
carried out in a manner which is essentially analytical as described 
above. This will not be done, however, because the modification 
is such that the resulting equations whose purpose is that of 
Equations (15) and (18) are too cumbersome to be of practical 
value. Thus the present discussion will be confined to a graphical 
solution of the problem. 


To evaluate /;, one notes that 
TAya7 
Ma Mr 
TAya1 i? Lat ~ Mr — Ma 4 Ma 
Ma “Mr 100 B, 


where M is the molecular weight of the indicated component. 
m, may then be evaluated in terms of b; and “aw, as 


100 





bh = 447 = 


yA =. 
as (21) 
100 — 74 AW. 


A useful substitution for i y4y in Equation (21) is 


_ Mw = 
aw Mw — Ma Ma 


Mu me Ct” «Pfs 





Equations (20) and (21) thus provide the means whereby the 
convenience of working on a weight % basis may be preserved 
until Equation (19) is required. 
It will be apparent that equations in terms of point f2 similar 
to those above for point i; may be stated. Thus 
f42 = meof¥2 + de. ; (23) 


fA2 _ Aaa be 


ss ., (24) 
100 — f4oan 


Mr 


Mr — Ma 
100 B, 


b, = 


As a transformation to mole % has occurred, it is appropriate 
to modify the approach of Equations (1) - (3) to 


Nis, + n; i4,y — ja 
N= = —.. (26) 
Nivy fA2 = 7 


where n; and n;,,, refer respectively to the number of moles of j 
added to iy. Unlike Equation (3), however, Equation (26) by 
itself is insufficient to locate i;y because Nj,,, 18 an unknown 
quantity. The solution to this difficulty may be found in the 
addition of another step to the procedure beyond that required 
for the use of Equation (3). This consists of the addition of a 
second working mixture of a known composition, such as h, to 
another sample of isy to achieve another point f;. Points fs 
and f, will not ordinarily lie on the same tie-line, except by 
coincidence. It will be apparent that equations in terms of f, 
similar to (23-25) above may be written. 


In the following, the restriction will be placed that neither 
j nor h is to contain any of that component whose concentration 
is used as the means of fixing the binodal points, i.e., in the 
present system, acetic acid. This restriction is necessary to 
obtain the resulting relationships in a usable form, as will be 
seen below. Comparable to Equation (26), 
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Ui, + Mh 


N’” = 
UNis yy 
where it is assumed that the sample of i,v to which h is added 
is V times as heavy as the one to which j was added. 
If to m4 moles of acid which are present in the mixture totalling 
n;,,, Moles there is added 7 containing 77 moles of toluene and 
Nw moles of water, then 


Let the amount of h added contain pn; moles of toluene and gnw 


moles of water. Then 


100 vn, 
fag Be ee eee eee (30) 
UNG, 1 + pnr + qnw 


and 


Ses (31) 


Ui 
Finally from Equations (34) and (37) 


Nn" - 
N 


pur + anw 


petcaiica ga aurea setae 33 
unr + ny) (33) 


Point i;jy may thus be located by striking straight lines from 
points j and h, which lines are to pass respectively through fr: 
and f; and intersect at a point on the tie-line containing i. The 
procedure has an obvious resemblance to that described for the 
weight % basis, except that Equation (32) must also be satisfied. 


It appears that, due to the larger number of steps involved, 
the mole % method is not capable of the best results obtainable 
on the weight % basis. The former would be preferable to the 
latter, however, when the complication of parallel tie-lines arises 
on the weight % basis. 
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Nomenclature 


A acetic acid 

Bb tie line intercept, wieght % basis 

b = intercept of tie-line, mole % basis 

Cc = plait point 

D mutual solubility of water and toluene 

E = extrapolated mixture, hypothetical 

F = mixture formed by addition of J(K) to upper (lower) 
phase of /, weight “% basis 
mixtures formed by addition of either 7 or h to upper 
phase of 7, mole % basis 
function (Equation (1)) of weights of J and upper phase 
of J 
function (Equation (16)) of weights of K and lower 
phase of / 
weight of a quantity of a component or mixture 
working mixture (contains no acid) used i in establishing 
composition of upper phase of 7, mole % basis 
mixture lying on tie-line containing binodal points which 
are to be established experimentally, weight % basis 
mixture lying on tie-line containing binodal points which 
are to be established experimentally, mole ©% basis 
working mixtures used in establishing composition of 
upper and lower phases of J, weight “% basis 
working mixture (contains no acid) used i in establishing 
composition of upper phase of 7, mole % basis 

5 lower phase of heterogeneous mixture 

M slope of tie-line, weight % basis 
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molecular weight 

slope of tie-line, mole % basis 

function (Equation (26)) of moles of 7 and upper phase 
of 7 

function (Equation (27)) of moles of # and upper phase 
of i 

moles of a quantity of a component or mixture 

weight of toluene (water) in quantity of h added to 
fraction of upper phase of i/weight of toluene (water) 
in quantity of 7 added to fraction of upper phase of 7 
reference mixture, prepared synthetically 

scanning mixture, prepared synthetically 

toluene 

upper phase of heterogeneous mixture 

weight of fraction of upper phase of i added to h/weight 
of fraction of upper phase of i added to j 

water 

working mixtures used in preparing R 

function (Equation (33)) of p, g, v and the moles of 
toluene and water in the quantity of 7 added to the 
upper phase of 7 
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NOTE TO THE EDITOR 


22 
A Gener alized Method for the Estimation m 
% : 18 
of Heat of Vaporization ns 
Te 
16 
CARL ST. PIERRE and CHI TIEN 14 
Department of Chemical Engineering, 
Essex College, 
Assumption University of Windsor, Windsor, Ont. 12 
a number of methods have been proposed for the prediction where a is the slope of the vapor pressure curve at the critical 10 
of heat of vaporization. However, most of these sugges- point. 
tions deal with a rather limited class of substances and in some Substituting (3) into (1), we obtain the expression for a at 
cases it requires at least one value of the heat of vaporization . the saturation condition: 
at a given — and os ete emai haa OWT:h = RZag — Zaz) Figure 
true sense. A recent paper by Kunte an raiswamy“ gave [aT;, + 0.0838 (a — 3.75) (36 + 67,7 — 427.) 
a correlation and claimed excellent accuracy. Unfortunately : : po 
there seems to have been some doubt about the validity of its Numerical values of Z,, and Z,, for specified values of ae see 
mathematical derivation and furthermore the computed value of and P, are given by Lyderson ct al® for values of Z, = 0.23, ceniian 


\. is found to be increasing with the increase of temperature 
(contrary to Trouton’s rule). 

The present communication aims at presenting a generalized 
method for estimating the heat of vaporization. As a starting 
point, the Clausius-Clapeyron equation as written in terms of 
compressibility factor is used. This is given as: 


h, RT,2(Zs¢ — Zsz) (AP, 
= _ = OT. oe . . eee (1) 


is P,, 
The heat of vaporization, \, refers to that at the saturation 
condition. Riedel’s expression for vapor pressure is used for 


Ap 
r 


0 
the evaluation of( ~), This is found to be: 


oT, 
ln P, — aln T, — 0.0838 (a — 3.75) 
(iT, — 35 ~ T+ 428 T,)....:€2) 
OP, 
and = 


aT, = P,{a/T, + 0.0838 (a — 3.75) 


(36/T,? + 67,5 — 42/T,)].... ..(3) 


0.25, 0.27 and 0.29. Values of a for different Z, can be found 
from Reference (3). With this information, values of (\/Tc), 
for a given Z, can be computed at various saturation conditions, 


For the values of (A/7T.) at temperatures other than the 
saturation conditions, the Trouton’s rule is used to account for 
the temperature effect, namely ; 


Xr 1 a r 0.38 
x [i 


From Equations (4) and (5), it can be seen that the numerical 
values of \/T, is a function of three parameters, 7,, P, and Z.. 
For any given value of Z., values of 4/7. are plotted against 
T, using P, as parameter. Figure 1 gives the result of calculation. 
As an indication of the reliability of this method, experimental 
values of \/7., for a number of compounds are compared with 
the values obtained from the graph. The accuracy is not as 


good as hoped. However, with the exception of chlorine, 


TABLE 1 


COMPARISON BETWEEN EXPERIMENTAL AND PREDICTED VALUES OF X 


Compound Ze 
| 


Argon) 87 
Chlorine® 238 
Carbon Dioxide® 28 194 
Hydrogen Chloride 27 188 


Hydrogen Sulfide 28 212. 


Carbon Tetrachloride™ 27 349 


(X) at boiling poi 


Boiling Point °K, 


Sulfur Trioxide” 26 316.. 


Sulfur Dioxide” 27 263 
Oxygen® 29 90 
Carbon Monoxide® 29 81 
Nitrogen 29 77 


Freon’ 27 242. 


Carbon Disulfide’ 29 319 
Hydrogen Cyanide 29 298 





Source of information: 


“Se UNmewN 


nt (cal./ml.) 


Experimental Predicted 


1,596 
4,878 
6,030 
3,860 
4,463 
7,170 
9,990 
5,950 
1,629 
1,444 
1,336 
4,850 
6,400 
6,027 


1,615 
5,660 
4,020 
4,530 
4,410 
7,400 
7,300 
5,880 
1,640 
1,385 
1,310 
5,100 
5,970 
7,350 


bro hd Ww 


NNO NN 











‘Selected Values of Chemical Thermodynamic Propeities, July 1953 
‘Chemical Engineers Handbook, 3 ed. McGraw-Hill Book Co, 1950 
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Nomenclature 
P, 


22.0 


reduced pressure 

gas law constant 

critical temperature 

reduced temperature 

compressibility factor 

slope of the vapor pressure curve at critical point 
heat of vaporization 


20.0 


Yea NSN 
vunnnnad 


18.0 . 
Subscripts: 

saturation state 
vapor state 
liquid state 


JR 
iow il 


16.0 
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Figure 1—Generalized chart for latent heat of vaporization. 
carbon dioxide, sulfur trioxide and hydrogen cyanide, it does 


give a fairly reasonable estimate. Table 1 summarized the 
results of comparison. 
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INDUSTRIAL SECTION | 


Pilot Plant Unit for the Continuous 
Methanolysis of Glyceride Oils’ 


P. M. JANGAARD? and F. A. VANDENHEUVEL? 


Efficient continuous methanolysis using an alka- 
line catalyst is achieved in the simple system described. 


At 100°C., in one hour, a volume of oil equal to 
six times the free volume of the reactor is converted 
to methyl esters. 


A three-step process for the production from herring oii of 
methyl ester fractions having specific properties was 
described by the authors in 1957"! The process comprised 
first, the production of mixed methyl esters by methanolysis of 
the oil; second, their fractionation by solvent segregation; and 
third, the fractional distillation of the segregated ester fractions. 
Since steps 2 and 3 lend themselves to continuous operation, it 
seemed desirable, both from the point of view of economy and 
of simplicity, to integrate the three steps in one continuous 
process. Continuous methanolysis was therefore investigated 
The reaction expressed by the following equation 

Ot 
Triglyceride + methyl alcohol 
had been known for many years as a typical example of alcohol 


glycerol + methyl esters 


ysis. As such it had been described and discussed in several 
articles which have been reviewed recently by Kaufmann et al’ 
and by Poré 


In many laboratories, the reaction was applicd to batch 
ty pe Operations only [he authors themselves had used ex 
tensively the conditions described by Bradshaw and Meuly“ 
whereby over 95% of a glvceride oil can be converted to 
methyl! esters in about one hour [he reaction was carried 
out at 60°C. using 75% exec methanol and 0.5% NaOH as 


the catalyst 


Ihe above conditions had been found very satisfactory in 
promoung the separation from the final reaction mixture of two 
ummuscible phases, the methyl ester layer and the glycerol layer, 
thereby unplifying the recovery of the product Increa ing the 
excess methanol beyond 100% was observed to interfere with 
this separation, and likewise, concentrations of the catalyst 


above 1% ere found to complicate subsequent operations 


Vas thought that a iniple olution to. the problem of 
j might consist of puliping mctered 


reavent through a heated channel allowing 
COCO eee ee Cees erecerseseeesseeeeeesseseceseseussecececscecsscssscecasscceseses 


Manuscript recetwed November 4, 1960; accepted March 20, 1961 
Fisheries Research Board of Canada, Technological Station, Halifax, N.S 
Department of Agriculture Animal Research Institute, Ottawa, Ont 
Contribution from the Fisheri« Research Board, Technological Station 
Halifax, N.S 
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Un reacteur simple 4 haut rendement horaire 
pour la méthanolyse continue en milieu alcalin est 
déerit. 


A 100°C. le volume d’huile convertie par heure 
égale 6 fois le volume libre du réacteur. 


sufficient time of contact for a satisfactory degree of methanolysis 
to develop in the mixture. 


Evidently, the efficiency of such a system would be increased 
by an increase of the reaction rate. For the reasons given above, 
means of attaining this result were restricted insofar as excess 
methanol and amount of catalyst were concerned. The effects 
of an elevation of the temperature, on the other hand, were 
unknown. 


The results of kinetic study of the reaction, which will be 


presented in detail elsewhere, indicate the following useful facts. 


\n appreciable increase in reaction rate can be achieved by 
operating at higher temperature. In Figure | is shown a curve 
expressing the relationship between the temperature and the 
tume required to promote 90% release of the glycerol from 
herring oil using 0.5% NaOH and 70% excess methanol. It 
was found that the final equilibrium composition did not change 
within the range of temperature explored. Under the conditions 
just described, a maximum release of the glycerol corresponding 


to 95%, was obtained in less than 7 minutes at 120°C. 


The reaction was also found to involve the consecutive 


equilibria described by the following equations. 
OH 
MeOH 
OH 
Diglyceride t MeOH 


lriglvceride + diglyceride methyl ester 


monoglyceride methyl ester 
OH 
MeOH methyl ester 


VMonoglyceride t vlycerol t 


\s these equations indicate, percentage Conversion to methyl 


ester will be higher than percentage glycerol release 


It was also found that continuous agitation of the reaction 
muxture did not decreas¢ the tine required to reach final 
equilibrium. Indeed, no agitation is necessary after an initial 
thorough mixing of the reagents 


In designing the sunple unit described below, we gave the 


reactor dimensions cstumated to correspond to an output of 2 


Su 
25g. 
methyl esters per minute at 110°C., in order to match the 
capacity of the solvent segregation unit used in the following 


step of oul proce ; 
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Figure 1—Curve expressing the relationship between the 

temperature and the time required to promote 90% release 

of the glycerol from herring oil using 0.5% NaOH and 
70% excess methanol. 


Allen and Kline have patented a procedure for the continuous 
methanolysis of glyceride oils‘, a general purpose continuous 
reactor has been patented, and a continuous laboratory unit 
for the conversion of acidulated soapstock using acid cataly st 
was described by Eaves et al. 


Description of the Unit 


As shown in Figure 2, a solution of the catalyst in absolute 
methyl alcohol is contained in tank A. This solution, and the 
oil stored in tank B are pumped simultaneously through the 
system, respectively by pumps C and D. The latter are actuated 
by a common adjustable speed motor drive E. The liquids are 
mixed in F where they form an emulsion. F is a shallow circular 
box containing a magnetic stirrer actuated by a motor-driven 
external magnet. The bolted lid of this box is fitted with a 
thick glass window allowing visual inspection for proper adjust- 
ment of the stirring speed. The latter should be sufficient to 
emulsify the mixture and maintain the box full of emulsion. 
The reaction mixture reaches preheating coil G, a 50-ft. length 
of coiled 14-in. internal diameter stainless steel tubing. This is 
followed by H, the reactor proper, a 50-ft. length of coiled 
34-in. internal diameter stainless steel tubing filled with small 
glass Raschig rings, the function of which is to promote sr 
heat exchange through agitation. Both coils are immersed 1 
oil bath I, which is ‘equipped with electric immersion me 
stirrer, temperature indicator and bimetallic thermoregulator. 
The latter, connected to a relay and switch combination, allows 
+0.5°C. adjustment of the bath temperature. ‘The temperature 
of the reaction mixture on entering the reactor is read on the 
recorder connected to thermocouple R This temperature 1s 
about 3°C. below the bath temperature. 

The reaction mixture then passes through cooling coil K, 
30 ft. of 14-in. internal diameter stainless steel tubing fitted in 
water jacket L, the flow being adjusted to bring the product to 
about 25°C. Spring valve M is adjusted to a pressure appreciably 
above the boiling point of the reaction mixture, 1e., 20 Ib 
above the vapor pressure of pure methanol at the temperature 
of the experiment. ‘Through valve N, the reaction mixture 
reaches centrifuge O (Sharples motor-driven super centrifuge) 
where the methyl! ester phase (outlet P) and the glycerol phase 
(outlet Q) are separated 


Operational Data 


\fter adjusting the bath temperature and the operating pres 
sure (valve M), three-stage metering pumps ¢ and I) are 
adjusted separately (stroke length) to deliver the required 
amounts of oil and catalyst solution at the selected motor speed 


The Canadian Journal of Chemical Engineering, August, 1961 


Figure 2—Pilot plant unit for the continuous methanolysis 
of glyceride oils. 


The liquids can then be pumped in simultaneously. 


A 100-ml. sample of product is collected at intervals from 
3-way valve N in a container containing enough mineral acid 
to neutralize the catalyst. Free glycerol is then determined by 
a standard method on the chloroform extracted mixture’. 
Constancy of the glycerol value is observed after 30 min. 
Valve N can then be connected to O, and the glycerol and 
methyl ester phases collected separately. 


Further processing of the methyl ester phase consists of 
vacuum stripping of the methanol it contains, and separating 
the small amount of soap which then precipitates. The glycerol 
phase is processed for the production of pure anhydrous glycerol 
according to a procedure which will be described in a separate 
article. 


Our work with this unit has consisted so far in establishing 
its efficiency under conditions of normal production. The data 
thus far collected therefore do not involve a wide range of 
conditions. : 


Dried oils were used throughout these experiments since it 
was established that moisture will interfere very markedly with 
the reaction at elevated temperatures. 


Some typical results shown in Table | have been obtained 
with seal oi] having the following characteristics: F.F.A., 1.5%; 
1.V., 146; total glycerol, 10.02%. With this oil, the maximum 
value observed for glycerol release was 93% in a batch operation 


after 24 hr. at 60°C. (100% excess methanol, 0.5°% NaOH) 


\ comparison of this result with those listed in Table 1 is 
legitimated by the simple rules governing equilibrium reactions 
It is well known that the final equilibrium composition 1s 
unaffected by catalyst concentration. Further, the equilibrium 
will be displaced by a larger excess of methanol towards higher 
levels of conversion Finally, it has been established, as stated 


TABLE 1 


GLYCEROL RELEASED FROM SEAL OW IN RELATION TO 
PEMPERATURE, TIME oF Contract, EXcEss METHANOL 
AND CATALYST CONCENTRATION 


contact MeOH 
niin 


Fotal flow | [ime ot} Excess | nioy Free glycerol 


ml. /min 


( 





above, that the degree of conversion was unaffected by the 
temperature used within the range of temperature explored. 
For these reasons, the 93% conversion obtained after 24 hours 
at 60°C., undoubtedly an equilibrium value, and obtained with 
the largest excess methanol applied in this study, must correspond 
to the maximum degree of conversion attainable with this oil 
within the range of conditions used. 

In Table 1, excess methanol is calculated as percentage over 
theoretical amount estimated from total bound glycerol. The 
catalyst used is expressed as percentage of oil weight. The 
time of contact is determined by dividing the free volume of 
the reactor (3000 ml.) by the total volume of reaction mixture 
flowing per minute. 

It appears that a degree of conversion approaching the 
maximum attainable is achieved within a time of contact pre- 
dictable from the kinetic data. The efficiency appears to fall 
to some extent when the temperature rises above 100°C. At 
this temperature the yield in methy] esters is about 260 gm./min. 
or somewhat in excess of the expected output. The volume of 
oil which can be treated in 24 hr. is 144 times the free volume 
of the reactor, an indication of the remarkable efficiency of this 
simple system. 
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Safety of Materials of Construction Used 


in Explosives Manufacture’ 


F. VAN ZEGGEREN? 


Hazards of materials of construction in machinery 
used for manufacture of explosives depend on the 
explosive which is being handled by the machinery. 
Thus one has to consider the combination material/ 
explosive as well as material and explosive them- 
selves, in assessing hazards. Several types of hazard 
ean occur viz., due to sparks, friction, impact or heat. 
For a given material the most likely factors that will 
influence its hazard depend on the type of hazard, 
e.g. for spark hazards important factors are: static 
pick-up of material and explosive, conductive proper- 
ties and spark confinements; for impact hazards: 
hardness, surface roughness, elasticity, presence of 
cavities, confinement offered and area of contact; for 
frictional hazards: hardness, friction coefficient and 
thermal conductivity. For purely thermal hazards, 
only flammability need be considered. 


In impact and friction the chemical action between 
material and explosive may play an important role. 
In all cases the sensitivity of the explosives themselves 
to the particular type of action, has to be taken into 
consideration. In accidental mishaps the type of haz- 
ard is often a combination of impact and friction, 
and requires a special test method. 


In this paper the literature on some aspects of 
all of the above hazards is reviewed. Experimental 
results of frictional and frictional impact tests are 
discussed from a theoretical point of view. A number 
of results of particular interest, e.g. for aluminum 
alloys, coated steel and non-metallic materials are 
presented. 


I" factories where explosive materials are being handled with 
machinery, the hazard associated with the contacting of 
explosive with the various parts of the equipment, will have to 
be extremely slight. Such hazards will depend on a variety of 
factors, viz., on the type of material in contact with the explosiv e, 
on the ‘explosive, and on a number of other conditions imposed 
upon the system by the type of hazard. Thus, unsafe situations 
can occur for different types of hazard, viz., due to sparks, 
friction, impact or heat. The specific safety requirements for the 
material to be used will, therefore, also vary. 

‘The mechanism operative in the initiation of an explosive is 
not yet fully understood and many conflicting data have been 


‘Manuscript received January 21; accepted May 380, 1961. 
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reported in the literature. Although in the last decade theoretical 
progress has been made, the choice of materials for machinery 
has still been largely empirical. 

A large number of test methods for materials have been 
devised, whose results have been interpreted as figures indicating 
degrees of hazard for those materials when used in a plant, but 
safety predictions are still rather uncertain, and progress along 
these lines would be welcomed. 

In this paper a survey will be given of the available literature 
on hazards of materials in contact with explosives and several 
test methods and apparatus will be discussed. A number of 
interesting data obtained for frictional and frictional impact 
hazards will be viewed in the light of recent theoretical advances 
in the field of friction safety predictions. 


Spark Hazards 


If an electrically conductive body is insulated from ground, 
it can act as a condenser or reservoir for the accumulation of 
electric charge. The close approach of another conductor may 
bring about ‘discharge in the form of a spark. An explosive 
material in the path of the spark may be ignited if it is sensitive 
enough and if the spark energy is high enough. Although this 
hazard cannot be completely excluded in explosive manufacturing 
machinery, the danger can be decreased considerably by adequate 
grounding of the conductive materials. For non-conductive 
materials the problem is magnified many times, since grounding 
is not a sufficient cure. 


Also serious is the generation of static on the explosive 
itself. Continued rubbing leads to formation of a double layer 
by electron flow from the surface levels of one substance in 
contact with another. Rubbing of particles against the container 
wall, particularly if this wall is smooth, leads to very appreciable 
static on the particle surface. The trouble can be overcome 
by roughening the vessel wall or by coating it with particles of 
the explosive material to be handled. 


In all cases the sensitivity of the explosive itself to sparks 
is of paramount importance. Minimum spark energies have 
been determined for a large number of explosives®:®. It was 
found that the sensitivity increases with confinement of the 
explosive. For instance, unconfined nitroglycerine could not be 
ignited with a 12.5 Joule energy spark but when heavily 
confined, a 7000 ergs spark could ignite the explosive. Thus 
confinement of spark and/or explosive should be avoided. 


Prevention of static accumulation is the best safety pre- 
caution, This can be effected by introducing radioactive 
elements, by roughening contact surfaces”, or by coating 
the explosive particles with antistatic agents. The latter 
method is by far preferred if the coatings do not change the 
required properties of the explosive. 
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Impact Hazards 

When a weight is falling down on an explosive compound, 
supported by a block, called an anvil, there is a chance that 
explosion occurs. The sensitivity of an explosive to impact is 
measured by this chance, usually with a so-called drop test 
apparatus, and is expressed in — sepened to give a certain 
percentage of explosions (e.g., 0, 50 or 100%). ‘Taylor and 
Weale®:® found that the hardness of the anvil was of great 
importance, and that a softer material required a larger fall 
energy, and hence is safer to use. Lack of rigidity of the anvil 
causes the required energy to be high, e.g. W hen the anvil is 
supported on a rubber pedestal. Thus it is preferred to keep the 
different machine parts as flexible as possible. 

The composition of the alloy of which the falling weight 
consists, affects the hazard. Wohler and W enzelberg“® found 
that when they used aluminum foil to cover up their drop 
testing samples, the impact hazard was increased, due to react- 
ivity of the aluminum and the large heat of its reaction with 
explosive oxygen. 

If the explosive, when subjected to impact, is not confined, 
the free falling weight will cause viscous flow in a horizontal 
direction. Then friction between weight and explosive takes 
place and surface roughness of the weight will become of 
importance, a rougher surface giving the more hazardous condi- 
tion (see Friction Hazards). 

For the same energy, a larger area of contact between falling 
weight and explosive reduces the probability for ignition. 
Thus, in machinery where impacting occurs, the total energy 
should be distributed over as large an area as possible. 


Bowden and co-workers“! found that the presence of bubbles 
in explosives greatly increases their ignition probability. This 
was explained by a mechanism involving adiabatic compression 
of the bubble, and subsequent heating of the bubble and of the 
surrounding explosive to a temperature well above the ignition 
point. This mechanism also explains the behavior of a cavity 
in the striker, which causes the required energy to be much 
smaller than for a flat-ended weight. Thus cavities, and in 
general roughness, of materials of construction should be avoided. 

Relative sensitivity measurements have been made for a 
large number of explosives in contact with several materials, 
and for a large number of materials in contact with several 
explosives"”. Quantitative interpretations of these results are 
still completely lacking. 


Friction Hazards 

When two solids are rubbed together, the frictional energy 
is dissipated mainly in the form of heat and is concentrated in 
the regions of contact, which are the summits of the highest 
irregularities of the respective surfaces“*. Thus hot spots of 
the order of 1000°C. can be formed. These high temperature 
flashes last for times of the order of a millisecond. When temp- 
erature and duration are high and long enough these hot spots 
can cause explosion in a lay er of explosive between the two 
rubbing surfaces. Experiments on this aspect have been carried 
out by Bowden, Stone and Tudor“ with nitroglycerine. The 
temperature rise in rubbing is given by: 


= uWo 
~ 4a J (ki + ke) 


where yu is the coefficient of friction, W the load between the 
surfaces, v the velocity of sliding, a the radius of the circular 
region of contact, J the mechanical heat equivalent and k, and ke 
the thermal conductivities of the rubbing solids. The area of 
contact is determined by the yield pressure (p,) of the softer 
of the two materials, and since the value of the hardness (H) 
is usually a few percent less than the corresponding value of Pon 
the above expression becomes: 


uv 


AT ~ 
ki + ke 
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If one of the two materials has a much larger thermal conductivity 
than the other one, we get 


VH 
AT ~ Mw 
ky 
Thus, one might assume that there exists a correlation between 
the property uVH/k;, of a material and its relative hazard when 
in frictional contact with an explosive. 


Dempster, arguing along the lines indicated above, showed 
that this was indeed the case for materials such as gunmetal, 
brass and steel, but aluminum gave too high a hazard value®), 
Dempster carried out his experiments with the help of a friction 
wheel. Gurton, Evans and Boyd carried out similar tests, for 
a large number of alloys, and when one examines their results, 
expressed i in minimum speed of wheel turning to cause ignition 
of nitroglycerine, one finds a very good correlation (0.9( 0) with 
the proposed relation, again with three aluminum alloys as 
exceptions“, More recently, Dempster has tested an additional 
50 alloys in the same manner, and rated them in order of increas- 
ing hazard. The exceptional behavior of aluminum in most cases 
has been attributed to the chemical interaction between aluminum 
and the explosive’s oxygen, thus giving a hotter “hot spot’’ and 
consequently faster expansion of this hot spot and greater ease 
of initiation. 


At our laboratory, we carried out a number of hazard tests 
for purely frictional hazard, for 18 metallic and non-metallic 
materials. The friction wheel used (Figure 1) consists of a 
10-12 in. disk of a material A to be chosen, which is turned by 
a two h.p., d.c. motor, the speed of which can be controlled by 
a Variac speed control box, whose settings correspond with 
speeds from 0-1700 r.p.m. The motor is bolted onto a rigid 
steel frame and a base which rests on hard rubber pads in order 
to keep vibration to a minimum. On the disk rests a rounded 
\4-in. diameter pin of another material B, 4-5 in. from the centre 
of the shaft. The load on the pin can be varied from one ounce 
to five pounds. Diametrically opposite the pin a fine paint 
brush applies a few drops of the liquid explosive under investiga- 
tion to the surface. On turning of the wheel the drops of explosive 
are smoothed out in a thin circular ring, which passes between the 
pin and the disk. The speed of the motor is increased slowly, 
until explosion occurs. Then the velocity of turning of the 
wheel is determined with a tachometer. This velocity (S) is 
indicative of the relative safety of the particular combination 
material A /explosive/material B chosen, a large velocity indicat- 
ing a smaller degree of hazard. 


Some of the results are presented in Table 1. The correlation 


je ; u Vi 
between minimum velocity and — was, however, zero. 
. ki + ke 


VARIABLE LOAD 
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Figure 1—Schematic diagram of the friction wheel. 
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TABLE 1 


A COMPARISON OF PHYSICAL PROPERTIES AND HAZARDS 
oF Vartous MATERIALS AS DETERMINED 
BY FRICTION WHEEL TESTS 


A. Combination: Material (pin) — Stainless Steel (disk) 
| 











u VH oer 

Material ;. + és S(cm./sec.) 
Glass | 9.3 2.9 
Stainless steel 8.3 6.5 
Mild steel 7.6 4.7 
Brass S.3 4.1 
Al-57SF a1 Ree 
Al-54SB 2:3 52 
Al-65ST-6 Eid 6.0 
Al-1SF £.3 5.4 
Al-2S1/2H 2 6.0 

B. Combination: Brass (pin) Material (disk) 
] 7 

Material uw a 

ki + ke (cm. /sec. ) 
Cupro-nickel 5.9 8.5 
Stainless steel $:3 72 
Brass I 4.9 15.5 
Al-57S1/2H 4.8 6.2 
A]-57SF 4.8 8.4 
Al-65ST-6 4.3 6.0 
Brass III 4.3 12.6 
Al-24ST-3 Sia 47 
Al-54SB 1.9 7.4 
Al-1SF 7 6.0 
Bakelite 6.5 >32 
Maple 2.0 >32 
Fir ye >32 


The explanation of this behavior probably lies in the false 
assumption of constancy of the properties u, H and k as friction 
is progressing. 
the disk occurs, and wu, H and k will change. However, repeated 
runs over the same area of contact, still yielded reproducible 
minimum velocity values. 

In view of these experiments, a pendulum apparatus such 
as described by Bowden and Yoffe” is probably preferable for 
measurements of pure frictional hazards, since there the temper- 
ature increase is almost instantaneous, and complications such 
as bulk heating and abrasion would not arise. The pendulum 
apparatus might be approximated by using a friction wheel with 
short tests of one or two seconds duration. 

In accidental mishaps, such as when a hammer drops on an 
explosive soiled floor, the hazard can best be tested by means of 
an apparatus simulating glancing blows. The torpedo frictional 
impact tester has found world wide application for this prupose. 
The apparatus, illustrated in Figure 2, consists of a brass 
V-shaped trough along an 8.7 Ib. steel torpedo with a nose 
hardened to 700 Brinell, can slide down a small anvil made of 
the material under investigation. The small square of material 
is clamped into a steel base plate which rests on a concrete pad. 
The trough is inclined under a variable angle, usually 45°. A 
smail amount of explosive is put on the surface of the awil and 
covered with a standard type sandpaper. The torpedo is then 
caused to slide down on to the sample from varying fall heights, 
and the 50% explosion point is determined (/,). 

Some of the test results are presented as semi-logarithmic 
plots in Figures 3, 4 and 5. The values of }, plotted in Figures 
3, + and 5 are quite reproducible. Usually, values are within 
0.2 inches from the average value given, and the maximum 
deviation obtained is +0.5 inches. 

One observes that, within the group of aluminum alloys, 
the relation uWH/k correlates very well with the vertical fall 
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As the pin slides over the disk, abrasion of 
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Figure 2—Schematic diagram of the torpedo frictional 
impact t tester. 


height value, for PETN (omnes ery ‘dialed tetranitrate), R- NG 
(regular —— cerine) and for 75G-NG, which is a nitrated 
mixture of 75% ethylene glycol and 25% glycerol. All three 
explosives ws were commercial grades. No attempts were 
made to purify the materials, since it was thought desirable to 
test the explosives that are actually being used in an explosives 
plant. 
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Figure 3—Comparison of physical properties and hazards 
of various materials in contact with PETN, as determined 
by torpedo friction tests. The circles denote aluminum 
alloys, viz.: 
@ ISF; © S50SF; ® 2S-1/2H; © 57S-1/2H; @ B54S; 

@ 24ST-3; @ 65ST-6; @57SF; @ 17ST-4. 
The squares denote other metals, viz.: 
@ Brass III; @ Cupro-nickel; BM) Brass 11; Nickel; 

Brass I; (W Steel. 





Figure 4—Comparison of physical properties and hazards 
of various materials in contact with R-NG, as determined 
by torpedo friction tests. The symbols denote: 


A Al 2S-1/2H; A A150 SF; & Al 65 ST-6; & Al 17 ST-4; 

® Brass III; @Brass I1; () Brass 1MSteel and Nickel; 

@ Cupro-nickel; @ Teflon; © Maple; @ Oak; © Maso. 

nite; @ Painted steel; © Transite; @ Bakelite; ® Nylon; 
@ Hard rubber; @ Tufnol. 


Within the group of other metals, the correlation is poor in 
Figure 3, and fair in Figures 4 and 5. In Figure 3 a dashed line 
instead of a solid one is used to express the uncertainty. It is 
interesting to notice the reversed positioning of cupro-nickel 
and brass I in Figure 3, and in Figures 4 and 5. When PETN 
is used cupro-nickel is the safest of the non-aluminum alloys 
and brass I is the worst. However, with R-NG and 75G-NG, 
brass I is the safest and cupro-nickel almost the worst. This 
type of behavior was also noticed by Corbett, of African 
Explosives and Chemical Industries*®. One observes from 
Figures 3, 4 and 5 that, although aluminum alloys are more 
hazardous than expected from their function uvH/k, several of 
them are still safer to use than, say, brass, e.g. Al 1 SF. 

For non-metallic materials (see Figures 4 and 5) the correla- 
tion is poor, mainly because of two materials: painted steel and 
masonite. In the case of painted steel, the exceptionally low 
value of the vertical fall height indicates a large hazard, larger 
even than for ordinary steel. It is thought that here the paint 
film is broken and that the hardness of the steel itself, combined 
with the low thermal conductance of the paint creates the high 
temperature which is thought to be responsible for the hazard. 
The behavior of the masonite stands still unexplained. 

Thus, although some of the test results can be explained 
with the theory based on the influence of physical material 
properties, predictions still appear to be largely unsatisfactory. 


Thermal Hazards 


In this category one need only consider the flammability of 
materials of construction, e.g., when several parts of the machin- 
ery are subject to continued friction, and one part inflames, the 
flame may propagate to parts in contact with the explosive. 
From this point of view, wood, although an ideal material to 
use from the point of view of sparks, friction or impact (see 
above), constitutes a real hazard. Fortunately, this hazardous 
situation arises from a very slow process, which can usually be 
detected in reasonable time by the machine operator. 
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Figure 5—Comparison of physical properties and hazards 

of various materials in contact with 75G-NG, as determined 

by torpedo friction tests. The same symbols are used as in 
Figure 4. 


Nomenclature 
AT 


temperature increment caused by friction, degrees C. 
coefficient of kinetic friction, dimensionless. 

load between two surfaces, dynes. 

velocity of sliding of two materials over each other, 
cm./sec. 

radius of circular regions of contact between two mate- 
rials, cm. 

mechanical heat equivalent: 1 cal. = 4.18 X 107 ergs, 
thermal conductivity of material, cal./cm. sec. degrees C, 
hardness of material, Brinell hardness units. 

yield pressure of material, kg./mm.? 

linear velocity of rotation of friction wheel, cm./sec. 
vertical fall height in torpedo tests, inches. 
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